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ABSTRACT 


We  consider  the  problem  of  reconstructing  the  flexural  rigidity  r(x) 

and  the  density  p(x)  of  a  beam.  The  unknown  beam  is  assumed  to  have  a  free 

left  end  and  a  clamped  right  one.  The  data  consist  of  the  displacement  and 

angle  of  the  center  line  of  the  free  left  end  following  an  initial  impulse. 

The  information  content  of  this  seismogram-like  impulse  response  is 

equivalent  to  three  spectra  { cu  },  {v  },  {y  }  and  two  gross  constants  F-, 

n  n  n  i 

?2‘  This  data  do  not  specify  the  structure  of  the  vibrating  beam  uniquely, 
but  rather  a  class  of  beams.  All  the  beams  in  this  class  share  the  same 
structure  over  that  portion  of  their  length  which  is  actually  set  in 
motion;  they  can  differ  over  the  portion  which  is  stationary.  A  method 
for  constructing  r(x)  and  p(x)  is  presented.  It  consists  of  two  steps: 
first  p(x)  and  r(x)  are  determined  over  a  small  interval  (0,x)  adjacent 


to  the  free  left  end.  Next,  this  known  portion  of  the  beam  is  stripped- 
off  and  the  response  of  the  resulting  truncated  beam  is  computed  via  the 
initial  data.  The  procedure  is  then  repeated.  Finally,  the  question  of 
the  existence  of  a  solution  is  discussed.  More  specifically,  conditions  on 
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INTRODUCTION 


This  paper  is  devoted  to  a  review  of  the  current  status  of  the 
inverse  problem  for  a  vibrating  beam.  This  inverse  problem  consists  in 
reconstructing  the  flexural  rigidity  r(x)  and  the  density  p(x)  as  a  function 
of  position  x  from  data  associated  with  natural  frequencies  of  vibration  of 
the  beam. 

The  reconstruction  of  r(x)  and  p(x)  is  made  without  assuming  any 
a-priori  knowledge  about  the  structure  of  the  beam.  Consequently,  this 
approach  differs  from  the  traditional  approach  for  the  geophysical  problem 
dealing  with  the  reconstruction  of  the  internal  structure  of  the  Earth  from 
data  on  toroidal  and  spheroidal  oscillations.  There,  the  spectral  data  are 
used  to  oorrecz  a  model  which  incorporates  knowledge  from  travel-time, 
surface  waves,  etc. . .  Thus,  the  Backus-Gilbert  technique,  which  is  ideally 
suited  for  such  correcting  tasks,  is  not  suitable  here. 

The  above  mentioned  geophysical  problem  is  a  primary  motivation  for 
the  investigation  of  the  inverse  problem  for  a  vibrating  beam.  Indeed,  if 
one  neglects  the  gravitational  force,  the  rotation,  the  oblateness,  etc... 
and  considered  the  Earth  as  an  elastic,  radially  stratified  sphere,  then 
its  normal  modes  are  of  the  two  types  previously  alluded,  namely  the 
spheroidal  modes  and  the  toroidal  modes.  The  latter  are  governed  by  a  2nd 
order  equation  of  Sturm-Liouville  type  whereas  the  former  are  governed  by 
a  4th  order  system.  The  bulk  and  shear  moduli!  as  well  as  the  density, 
which  characterize  the  elastic  properties  of  this  idealized  Earth  are  assumed 
unknown.  The  inverse  problem  consists  in  retrieving  these  characteristics. 

In  order  to  understand  this  difficult  inverse  problem,  it  might  be  helpful 
to  consider  simpler  inverse  problems,  and  in  particular,  simpler  inverse 
problems  associated  with  4th  order  equations  since  the  inverse  Sturm-Liouville 


problem  is  well  understood.  The  inverse  problem  for  the  vibrating  beam  is 
an  ideal  candidate. 

This  problem  can  also  be  viewed  as  an  example  of  an  inverse  problem 
associated  with  a  class  of  differential  operators  of  order  2n,  namely 

r  ,  , .  n  d  d  d  d  d  d 

2n  o  dx  1  dx  dx  n  dx  dx  1  dx  o 

where  all  the  a's  are  positive.  For  certain  boundary  conditions,  the  Green's 
functions  associated  with  these  operators  are  oscillating  kernels.  The 
theory  of  such  kernels  plays  an  important  role  in  the  solution  of  the  inverse 
problem. 

The  outline  of  the  paper  is  as  follows:  in  §1,  we  shall  introduce 
the  terminology  associated  with  the  five  eigenvalue  problems  which  will  be 
required  in  the  sequel.  Various  properties  of  fundamental  solutions  will 
also  be  given.  The  question  of  the  uniqueness  of  the  inverse  problem  is 
taken  up  in  §2.  We  shall  show  that  the  impulse  response,  which  contains 
information  equivalent  to  three  spectra  and  two  gross  constants,  determines 
a  class  of  beams.  The  beams  in  this  class  have  the  same  massive,  flexible 
front  end  of  length  L,  and  a  weightless,  infinitely  rigid  back  end  of 
arbitrary  length.  By  selecting  that  non-pathological  beam  without  any 
back-end  we  shall  achieve  uniqueness.  The  actual  construction  of  this  beam 
is  presented  in  §3.  The  construction  procedure  requires  a  discretization 
of  the  vibrating  system:  this  is  accomplished  by  replacing  the  beam  by  a 
series  of  massive  clothespins  of  various  stiffness  connected  to  each  others 
by  weightless,  infinitely  rigid  rods.  Formulas  for  the  characteristics  of 
the  clothespin  nearest  to  the  excitation  are  given.  Once  this  clothespin 
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is  determined,  it  is  "stripped-off"  the  structure;  the  spectra  for  the 
reduced  beam  are  found  in  terms  of  the  given  data.  The  same  procedure  is 
repeated  until  all  of  the  characteristics  are  determined.  This  procedure 
fails  if  the  given  data  is  not  bona  fide  data.  The  question  of  what 
constitutes  bona  fide  data  is  also  discussed. 
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PRELIMINARY  RESULTS 


We  shall  be  concerned  with  various  solutions  of  the  equation 

Hx*  =  “2  pOO  y  »  0  <  x  <  L  .  (1.1 

The  flexural  rigidity  r(x)  and  the  density  p(x)  are  non-negative 
functions,  i.e., 

r(x)  >  0  , 
p (x)  >0  . 

In  the  section  dealing  with  the  uniqueness  of  the  inverse  eigenvalue  pro¬ 
blem,  we  shall  assume  that  r(x)  and  p(x)  are  differentiable  and  in 
fact  that  they  possess  four  derivatives.  These  requirements  are  very 
restrictive  and  can  be  relaxed  in  places. 


We  also  find  it  convenient  to  write  (1.1)  as  a  system  of  first  order 
differential  equations,  namely 


de  _  1 
dx  ~  r  T 

dt  _ 

dx  "  ‘  x 

=  -  a)2  p  y 

dx 


(1. 


The  variables  9  ,  t  and  x 
the  slope  of  the  central  line,  t 


have  a  simple  physical  meaning:  9  is 
is  the  stress  and  x  the  moment. 
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Several  fundamental  solutions  of  (1.1)  will  be  useful  in  the  sequel. 

In  particular,  we  shall  define  a>2)  and  u>2)  thus: 

$(0,  u>2)  -  1  =  4>'(0,  oj2}  =  ■$ '  1  (0,  oj2}  =  4>'  1 1  (0,  u) 2)  =  0  (1.3a) 

*(0,  oj2)  =  VtO,  w2)  -  1  =  ^'(0,  w2)  =  *,,,C0,  w2)  =  0  (1.3b) 

when  primes  denote  derivatives  with  respect  to  x  .  <j>(x,  oj2)  and  ij;(x,  w2) 

are  entire  functions  of  oj2  of  order  .*  To  prove  this  assertion  one 
would  show  (i)  that  the  Taylor  series  expansions  of  $  and  are  convergent 
throughout  the  entire  complex  w2-plane  and  (ii)  that  $  and  ^  are 
dominated  by  entire  functions  of  w2  of  order  ^  .  The  reader  is  referred 
to  Titchmarsh  (1962)  p.  6-11  where  a  similar  proof  is  given  for  the  Sturm- 
Lionville  case. 


see  Boas  (1954)  p.  8  for  the  definition  of  "order". 
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By  replacing  s  by  -  co2  ,  we  see  that  the  entire  functions  <*>  and  ^ 
are  indeed  of  order  ^  . 

$  and  ^  are  particularly  well  suited  for  studying  the  vibrations 
of  a  beam  with  a  free  left  end,  i.e.  a  beam  satisfying  the  boundary  con¬ 
ditions 

y1  1  =  Cry11)1  =0  at  x  =  0.  (1-6) 

Insofar  as  the  right  end  is  concerned,  we  shall  exclusively  deal  with  the 
clamped  case,  i.e. 

y  =  y 1  =  0  at  x  =  L  .  (1.7a) 

To  that  effect,  we  introduce  two  additional  fundamental  solutions  of  (1.1), 
say,  u  and  v  such  that 

u(l,  a,2)  =  u 1  (L ,  (o2 )  =  u"(L,  <o2)  -  1  =  u'"(L,  <o2)  =  0  ,  (1.7b) 

v(L,  co2)  =  v'(L,  co2)  =  v 1  1  (L ,  co2)  =  v"'(L,  <o2)  -  1  =  0  .  (1.7c) 

Needless  to  say,  u  and  v  are  also  entire  functions  of  co2  of  order  1  . 

The  Wronskian  W(x,  co2)  which  is: 

41  u  v 

u'  v' 

0"  u ' '  v' 


(1-8) 
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satisfies  the  following  differential  equation: 


^  [r2(x)  W(x,  u.2)]  =  0  .  (1.9) 

Consequently 

r2(x)  W(x,  u2)  =  r2(0)  W(0,  to2)  ,  (1.10a) 

=  r2(L)  W(L,  u)2)  .  (1.10b) 

But 

W(0,  w2)  =  u"  (0,  oj2)  v' 1 '  (0,  u,2)  -  u' 1  ‘  (0,  «2)  v"  (0,  w2)  (1.10c) 

and 

W(L,  to2)  =  ( L ,  q'(L,  uj2)  -  o'(L,  u^)  -j;(L,  u2)  .  (l.lOd) 


Therefore  the  Wronskian  W(x,  u2)  vanishes  whenever  m  coincides  with  the 
eigenfrequency  oon  of  the  beam  in  the  free/clamped  configuration. 

The  eigenvalue  problem 


K' 

)"  = 

un  P  ^n 

y*  - 

K' 

-  0  , 

for  x  =  0  , 

►  (1.11) 

= 

=  0  , 

for  x  =  L  ,  j 

will  be  the  central  eigenvalue  problem  of  the  paper.  However,  other  eigen¬ 
value  problems  will  occasionally  be  needed.  They  differ  from  (1.11)  only 
insofar  as  the  left  boundary  conditions  are  concerned:  the  right  end  is 
always  as  in  (1.11),  i.e.  clamped.  For  instance,  if  the  left  end  were 
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clamped,  then  we  would  have: 


The  eigenvalue  problems  (1.11),  (1.12),  (1.13)  and  (1.14)  are 
self-adjoint  and  physically  meaningful.  The  next  pair  of  eigenvalue  pro¬ 
blems  are  neither.  However,  they  do  enter  in  our  analysis: 


ton 

'!"  ■ 

d  y 
Jr, 

i  ’ 

' 

M 

c 

to(' 

)•  ■ 

--  0 

at 

O 

II 

X 

! 

>  (i. 

,15a) 

*n  = 

y' 

7n 

■  0 

at 

x  =  L  ; 

j 

K 

’)••  = 

Si 

o  y 
•'n 

i  » 

i 

yh  = 

y"  = 
Jrt 

0 

at 

X  « 

0  , 

\ 

(i. 

,15b) 

yn  = 

y'  = 

Jn 

0 

at 

X  = 

L  . 

i 

j 

(1:15b)  is  the  adjoint  of  (1.15a):  hence  they  have  the  same  eigenfrequencies. 

More  general  boundary  conditions  could  have  been  considered,  but  for 
the  sake  of  presentation,  I  have  restricted  myself  to  the  above  simple  con¬ 
figurations.  The  Table  summarizes  the  notations  and  is  convenient  for  future 
reference 


Left  Boundary 
Conditional 

Eigenfrequencies 

. 

Configuration 

y  =  y*  =0 

Xn 

clamped 

. 

y  =  y"  =0 

% 

supported 

y  =  (ry1 1 )•  =  0 

vn 

non  self-adjoint 

y>  =  y"  =0 

%  ' 

non  self-adjoint 

an 

Rayleigh 

y"  =  (ry")1  =  0 

wn 

free 

_ 1 
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Let  us  return  to  our  fundamental  solutions  <t>,  ■?,  u  and  v  and  record 
their  values  for  m2  =  0  : 

<p(x,  0)  =  1  , 

4>(x,  0}  =  x  , 


(1.16a) 

fL  *>■ 

(1.16b) 

lL  nh  • 

(1.16c) 

fL  tx-x;ixcT~u  *•  ■ 

(1.1 6d ) 

As  a  result,  we  deduce  from  (l.lOd)  that 

W(L ,  0)  *  1  O-V) 

Making  use  of  Hadamard's  Factorization  theorem  for  entire  functions 
(Boas  1954,  p.  22)  for  W(L,  w2)  which  is  also  an  entire  of  order  j  ,  we 

can  write 

W(L,  u,2)  =  W(L,  0)  n  0  -  fr  )  •  (1J8a) 

m=l  m 


or,  on  account  of  (1.17), 

W(L,  a)2 )  =  ”  (1  -  )  • 

m=l  “m 


(1.18b) 


The  asymptotic  behavior  of  the  eigenfrequencies  {u>nl  for  n  large, 
which  we  shall  recall  presently,  guarantees  the  convergence  of  the  infinite 
product.  Inserting  (1.18a)  in  (1.10b),  we  get 
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(1.19) 


W(x,  ID2)  = 


n  (1  ■  ~~r  )  • 


The  fundamental  solutions  u  and  v  will  not  be  used  in  their 
present  form.  Rather  they  will  enable  us  to  define  two  new  fundamental 
solutions,  viz 


and 


<t(x,  w2)  =  r(x) 


0  U  V 

4,'  u1  v' 

ij>'  1  u' 1  v 1 1 


T (x ,  W2)  =  r(x) 


<6  u  v 

4>  *  u’  v’ 

4> '  *  u‘ 1  v' ' 


(1.20) 


(1.21) 


It  is  a  simple  matter  to  check  that  $  and  ?  are  indeed  solutions  of  (1) 
as  well  as  that 


<t>''(0,  id2)  =  $(L,  id2)  =  <5>’(L,  id2)  =  0  ,  (1.22) 

(H'"),(0,  id2)  =  ?(L,  id2)  =  (L,  id2)  =  0  .  (1.23) 


We  give  the  asymptotic  behavior  of  #  and  <?  obtained  once  again  by  means 
of  a  WKBJ  approach 


4>(x, 


-s)  % 


1 


(1.24a) 


•  {[cos  Tjr  sinh  w  +  sin  w  cosh  w]  s i n ( X  -  w)  sinh(X  -  nr] 

cos -o  cosh  or  [sin(X  -  w)  cosh(X  >  «)  -  cos(X  -  «■)  sinh(X  -  -w)]} 
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H'(X,  -S)  % 


(1.24b) 


•  {2  cos  ©  cosh  ©  sin(X  -  w)  sinh(X  -  ©  ) 


+  [cos  ©  sinh  ©  -  sin  ©  cosh  ©][sin(X  -  ©)  cosh(X  -  ©)  -  cos(X  -  ©)  sinh(X  -  ©)]} 


where 


© 


dx  , 


and  X  is  given  in  (1.5).  Finally,  by  substituting  (1.4a)  and  (1.4b)  in 
(l.lOd)  we  can  deduce  that 


r2(0)  W(0,  -s)  -v,  •  r(0)  r3(L)l  2  (cos2  ©  +  cosh2  ©}  (1.25) 


In  view  of  (1.22),  we  can  use  the  function  $  to  solve  the  eigenvalue 
problems  (1.11),  (1.13)  and  (1.15b).  In  particular,  from  (1.13)  we  deduce 
that 

$(o,  u)2 )  =  <f(o,  o)  n  (l  -  -V  ) 

m  ym 

and,  making  use  of  the  formulas  (1.16)  which  express  the  behavior  near 
x  =  0  ,  we  see  that 

<t(0,  uj2)  =  r2(L)  Fj  n  (1  -  *2  )  (1.26) 

m  vm 
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0.27) 


(1.28) 


0.29) 


Finally,  by  similar  means  we  can  show  that 

*'(0,  w2)  =  -r2(l)  Fj  n  (1  -  )  .  (1.30a) 

m  m 

r(0)  «'"(0,  cc2)  =  r2(L)  n  (1  -  4  )  •  (1.30b) 

m  “m 

The  functions  $  and  f  are  linearly  independent  as  long  as  w  differs 
from  an  eigenvalue  wn  of  the  beam  in  the  free/clamped  configuration.  In 
this  case  we  have 

4>(x,  +w2)  =  kn  v(x,  u>2)  .  (1.31) 
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The  constants  of  proportionality  k.  ,  can  be  deduced  by  specializing 
(1.31)  to  x  =  0  and  using  (1.26)  and  (1.28): 


o 


(1.32) 


To  simplify  the  notations,  we  shall  define 

Tn(x)  =  f(x,  co2)  .  (1.33) 

Thus  {'t'  ,  to  }  are  the  eigensolutions  of  (1.11).  We  can  also  represent 
n  n 

the  eigenfunctions  of  (1.11)  in  terms  of  $  and  ip  .  Thus  we  have  another 
formula  akin  to  (1.31): 


ip(L,  gj2  ) 

*(x’  “n5  •  *(x’  wn}  =  K  ‘ 


(1.34) 


It  should  be  noted  that  from  the  definitions  (1.20)  and  (1.21)  of 
4>  and  f  we  can  deduce  that 


<t>''(L,  ui2)  =  r‘(L)  i>(L,  to2) 
T1 '  (L,  u>2}  =  r'  (L)  <p(L.  w2) 


As  a  result 


<p(L,  w*)  1  (L,  u2) 

<i> CL ,  go2)  t'1 '  (L,  oj2)  ^n 


(1.35) 


i 

i 


l 

t, 

ii 


and  so  (1.34)  can  be  written  as  follows: 


kn  <|>(x,  u£)  -  i|/Cx,  tuj)  =  k*  ^(x)  .  (1  •  34' ) 


In  order  to  find  the  values  of  k*  ,  we  set  x  =  0  and  use  (1.28): 


k; - = - ?!-  •  0-36) 

r2(L)  F,  S  (1  -  4  ) 
m  m 


We  consider  next  the  asymptotic  form  of  the  eigenfrequencies  of  the 


five  problems.  It  is  easy  to  see  that 

♦  (0,  Xj)  r{0,  xj)  -  *'(0,  X*)  y(0,  Xj)  =  0  (1.37a) 

♦(0,  uj)  =  0  (1.37b) 

f(0,  vj)  =  *'(0,  v2)  =  0  (1.37c) 

»(0,  cj)  =  0  (1.37d) 

(rV1)'  (0,  u.2)  =  f(0,  «2)  =  0  (1  •  37e ) 


As  a  matter  of  fact,  (1.37b)  and  (1,37c)  are  equivalent  to  (1.26),  (1.28) 
and  (1.30).  By  replacing  $  and  v  in  the  above  equations  by  their 
asymptotic  forms  as  given  in  (1.24a)  and  (1.24b)  we  deduce  that 

A*  t,  (n  +  H  )4  ttVs4  "I 

U2  *  («♦%)* 

^  n4  ir4/=4 

(n  -  H  )4  tt4/h4 
oj2  ^  (n  -  h  )4  tt4/e4 


(1.38) 


1.12 


where 


H  =  fL  C°-  )%  dx  (1.39) 

Jo 

Other  important  properties  of  the  eigenvalues  will  be  needed.  For 
instance,  we  shall  exploit  the  fact  that  the  eigenvalues  xn,  un,  vn,  an 
and  u)n  are  all  simple  eigenvalues.  This  result  is  a  consequence  of  the 
theory  of  integral  equations  with  oscillating  kernels  [Gantmacher  and  Krein, 
1950].  We  can  only  give  the  briefest  outline  here.  The  first  step  is  to 
convert  the  eigenvalue  problems  (1.11),  (1.12),  (1.13),  (1.14)  and  (1.15) 
into  integral  equations  of  the  form 

rL 

y(x)  -  u,2  G(x,  t)  P(t)  y(t)  dt  .  (1.40) 

JO 

If  the  kernels  G(x,  t)  ,  which  are  the  Green's  functions  for  the 
various  problems,  satisfy  the  following  properties: 


(i)  G(  J1  •*'  *n  )  e 

ti  ...  t 


G(xj,  tj)  ...  G(xj ,  tn) 
G(xp,  t  )  ...  G(xn,  tn) 


>  0 


(1.41) 


for  all  partitions 


Xj  <  x2  <  . • •  <  x 

0  <  <L,  n  =  1 ,  2 , 

tj  <  t2  <  •  •  •  <  t^ 


and 


Xi  ...  x 

(ID  G(  vn  )  >  0 


Xi  . . .  x. 


(1.42) 


for  all  partitions  and  all  values  of  n  ,  and 


(iii)  G(x,  t)  >  0 


for  0  <  x,  t  <  L 


(1.43) 


then  the  kernel  G(x,  t)  is  said  to  be  oscillating.  For  such  kernels, 
Gantmacher  (1936)  [see  also  Gantmacher  and  Krein,  1950]  showed  that  the 
eigenvalues  of  (1.40)  are  simple. 

Gantmacher  and  Krein  (1950)  have  shown  that  the  Green's  functions 
for  four  out  of  the  five  eigenvalue  problems  of  interest  to  us  are  indeed 
oscillating  kernels.  These  are  the  four  self-adjoint  problems  (1.11),  (1.12), 
(1.13)  and  (1.14):  thus  *n,  up,  ap  and  ajp  are  simple  eigenvalues. 

The  non-self  adjoint  problems  (1.15a)  has  not  been  considered.  How¬ 
ever,  by  means  of  a  theorem  of  Karlin  (1971),  we  can  state  that  the  Green's 
function  for  that  problem  satisfies  1.41,  i.e.  is  "totally  positive"  in 
Karlin's  parlance.  The  fact  that  tjie  kernel  of  the  integral  equation  is  a 
"totally  positive"  Green's  function  implies  that  (1.42)  is  also  satisfied: 
the  proof  of  this  assertion  can  be  found  in  Krein  and  Finkelstein  (1939). 

Thus  we  only  need  to  prove  that  (1.43)  holds. 

It  is  a  simple  matter  to  check  that  the  Green’s  function  under  con¬ 
sideration  denoted  by  N(x,  t)  is  given  by 


N(x,  t)  =  n(x,  t)  -  jl  |  dz 


z(z  -  t) 

t  “W 


dz 


(1.44) 


where 


(z  -  x)(z  -  t) 


g(x,  t) 


X  r(z) 


(L  (z  -  xHz  -  t) 


T[FT 


dz  ,  x  >  t  , 


dz  ,  x  <  t 


(1.45) 


Consequently, 


1.15 


where 


x<)(t)  = 


z  dz 

HU 


+ 


[L  z2  dz  1 

JtTHTj 


0.49) 


Clearly 


t  <  x0(t)  <  L 

and  so  N(x,  t)  is  respectively  concave  downward  and  upward  over 
(t,  x0(t)}  and  (x0(t),  L)  .  Note  that 

N(x0(t),  t)  =  fl(x0(t),  t)  -  [x0(t)  -  t]  |  Z  ''"rfzj^  dz 

x0(t) 


i  .e. 

N(x0(t),  t)  =  jL  dz  >  0 

x0(t) 


and  so 

N(x,  t)  >  0  for  0  <  x  <  x q ( t )  . 

See  Fig.  1 . 

We  obtain  the  desired  result,  by  noting  that  over  the  interval  (x0(t),  L), 
N(x,  t)  must  lie  above  its  tangents.  The  x-axis  being  one  such  tangent, 

N  must  therefore  be  positive  in  this  interval.  Consequently 

N(x,  t)  >  0  ,  for  0  <  x,  t  <  L  , 

i.e  N(x,  t)  is  an  oscillating  kernel  and  hence  the  eigenvalues  vn 
are  simple  eigenvalues. 


1.16 


Consequently, 


1.15 


IMPULSE  RESPONSE  AND  UNIQUENESS 
OF  THE  INVERSE  PROBLEM 


The  impulse  response:  Let  us  consider  the  following  artificial  problem 


3X 


r(x)  =  0<X<L 

t  >  0 


(2.1) 


wi  th 


(0>  t}  =  h  (r(x)  ) (0>  =  0  * 


y(L»  t)  =  —  y(L,  t) 
y(x,  0)  =  n(x) 


=  0 


(2.2) 


After  making  a  Laplace  transform  where 


y(x,  -s)  =  {  e"st  y(x ,  t)  dt  , 
o 


we  can  write  (2.1)  as  follows: 


Cry « ' ) ' '  +  p  s  y  =  d  n 


(2.3) 


with 


y"  =  (ry ' ' ) '  =  0  at  x  =  0 

y  =  y'  =0  at  x=L 


(2.4) 


2.1 


... 


In  order  to  solve  the  above  boundary  value  problem,  we  introduce  the 
appropriate  Green's  function 


(2.5) 


(2.6) 


If  we  were  to  replace  s  by  -w2  ,  then  we  would  be  back  to  our  beam  which 
is  here  excited  by  a  force  n(x)  applied  at  time  t  =  0  .  The  case 

i(x)  =  s(x)/p (0) 

would  then  correspond  to  a  point  force  applied  at  the  force  left  end 
time  t  =  0  .  By  impulse  response,  we  shall  denote  the  measurements  of 
the  displacement  and  slope  of  the  center  line  at  that  left  end,  i.e. 


y(0 ,  -s)  =  £(0,  0;  s) 
8(0,  -s)  =  (0,  0;  s) 


(2.7) 


Now,  making  use  of  the  functions  $,  ip,  t  and  v  introduced  previously, 
we  can  check  that 


n(x,  z;  -s)  = 


i(x,  -s)  Hz,  ~s ) -  (x ,  -s)  Hz,  -s) 
r2(0)  W(0,  -s) v  7 


$ (z >  -s)  <?(x,  -s ) -ifr (z ,  -s)  f(x,  -s) 
r2(0)  W(0,  -s) 


(2.8) 


2.2 


Consequently 


y(0.  -s) 


and 


e(0,  -s) 


1 

'  r^ToT 


or,  in  view  of  (1.19),  (1.26)  and  (1.28) 


y(0,  -s)  =  F2  “ 
m=l 


1  +  s/u* 
_  m 

1  +  s/ uiji 
m 


(2.9a) 


0(0,  -s)  =  -  Fj  r? 

m=l 


1  +  s/,4 


(2.9b) 


The  impulse  response  contains  the  following  information: 

(i)  the  first  and  second  moments  Fj  and  F2  of  the 
"flaccidity" 

Cii)  the  three  spectra  {oip},  Tvn }  and  {yn}  corresponding 
respectively  to  the  free,  non-self  adjoint  and  supported 
boundary  conditions  at  the  left  end. 

Is  the  information  contained  in  the  impulse  response  sufficient  to  char¬ 
acterize  a  beam  uniquely?  This  is  the  question  we  wish  to  discuss. 

Beams  with  identical  spectra:  Let  us  say  that  the  beam  characterized  by 
r(x),  p (x ) ,  of  length  L  ,  clamped  at  the  right  end,  has  spectra  {XnJ, 
{%}.  ^vn^’  *CTn^  and  ^“n^  ‘  There  are  otfier  distinct  beams  with  the 
same  spectra.  Indeed,  the  two  parameter  family  of  beams  characterized  by 


I 


(2.10) 


r(x)  =  «  r(§  )5*  x)  , 
p(x)  =  6  pC|)'4  x)  , 

have  the  same  spectra.  (The  parameters  a  and  e  are  positive.)  The  lengths 
of  these  beams  are 

t  -  (  |  )’*  L  .  (2.11) 

The  first  and  second  moments  of  the  flaccidity  of  these  beams  are: 


(2.12) 


If  we  were  to  impose  the  additional  requirement  that  these  two  moments 
match  those  of  the  beam  under  consideration,  then  we  would  obviously  see 
that  a  =  3  =  1  .  Hence,  the  role  of  Fj  and  F2  is  to  set  the  scale  of 
the  beam  which  is  being  reconstructed. 


Pathological  Beams  with  the  same  impulse  response:  It  is  possible  to  con¬ 
struct  two  distinct  beams  with  the  same  impulse  response? 

Consider  a  beam  p(x),  r(x)  of  length  L  ,  clamped  at  the  right  end. 
Construct  a  second  beam  with  the  following  characteristics: 

pCx)  0  <  x  <  L 

pCx)  =  (2.13a) 

p  L  <  x  <  t 


2.4 


r(x)  0  <  x  <  L 


(2.13b) 


r(x) 


L  <  x 


In  other  words,  the  second  beam  is  identical  to  the  first  one  over  the 
interval  CO,  L)  and  then  has  constant  properties  over  the  section 
CL,  t)  .  Abusing  our  notations,  let  us  denote  by  $(x,  w2),  ^(x,  w2) 

two  fundamental  solutions  of  0-1)  which  satisfy  the  boundary  conditions 
at  x  =  0  appropriate  for  either  one  of  the  three  configurations  entering 
in  the  impulse  response. 

The  solution  over  the  interval  CL,  £)  is  written  in  terms  of  functions 

cos[(  aai  )h  (X  -  £)]  -  cosh[(  }*  Cx  -  t)l 
r  r 

and 

since  asi  )h  (x  -  t)]  -  sinhCC  )h  Cx  - 1)] 
r  r 

which  satisfy  the  clamped  conditions  identically. 

Piecing  the  solutions  in  CO,  L)  and  (L,  t)  by  requiring  that 

[y]  =  Cy ' 3  =  Cry”]  =  C(r^‘  =  o  C2.14) 

where  [  ]  denotes  the  jump  from  L  +  0  to  L  -  0  ,  we  obtain  the 
determinantal  equations  yielding  the  various  spectra,  viz.. 


2.5 


4>(L,o)  ) 


2 

^(L,u)  ) 


cos  Z 
-  cosh  Z 


sin  Z 
-  sinh  Z 


4> '  (L,o)  )  'll '  (L,oT) 


(r$")(L,uT)  (rD  (L,w  ) 


(rD'a.u)  (rip")  ' 


^  w5  {sin  Z 


+  sinh  Z} 


-X- 

-r2p  2  u{cos  Z 
+  cosh  Z) 


r  -p^w^{sin  Z 
-sinh  Z} 


-  h,  i 

—  to  {cos  Z 


-  cosh  Z} 


r 

-r  p  ‘  cj{  sin  Z 


+  sinh  Z) 


i  _3  3 

-r  p4u>2  {cos  Z 
+  cosh  Z} 


where 


z  =  |  u)'5  a  -  l> 

r 


2.6 


Now,  we  consider  the  limit 


^  -  0  (2.16a) 

r 

but  such  that 

r  p  =  a  .  (2.1 6b) 

This  means,  that  in  the  limit  the  added  section  is  made  up  of  a  weightless 
rod  of  infinite  rigidity.  Then,  dividing  the  last  column  of  (2.15)  by 

1  3 

rJ  p4  and  taking  the  limit  mentioned  above,  we  note  that  (2.15)  becomes: 


♦ 

r4>" 

(rVT 


P  0  0 

<P'  0  0 

rip'  -2au>  -to 

3 

(r<p")'  0  -2w2 


0 


or  better  still 


P 

p' 


0  . 


But  this  is  the  generic  determinal  equation  for  the  smaller  beam!  In 
addition,  since  the  flaccidity  of  the  second  beam  is  zero  over  the 


additional  interval  (L,  t) 


=  Fi 

and 

h  =  f2 

Thus,  for  the  clamped  right  end  case,  given  any  beam  with  impulse  response 
UnK  (vn},  {«■>„} »  Ft  and  F2  ,  we  can  construct  a  class  of  beams  with  the 
same  impulse  response,  by  simply  tacking  onto  the  original  beam,  an  infinitely 
rigid  weightless  rod  of  arbitrary  length. 

This  result  is  easy  to  understand  physically:  the  section  (L,  t)  , 
because  of  its  very  nature,  is  never  set  in  motion  by  the  initial  excitation 
and  hence,  never  felt. 

The  question  remains:  are  there  other  types  of  less  pathological 
beams  which  have  the  same  impulse  response?  The  answer  is  no. 

Toward  the  Uniqueness  Theorem:  We  assume  that  two  beams  ,  r^  and 

(21  (21 

d '  ,  rv  '  have  the  same  impulse  response.  Since  the  information  about 

the  length  of  the  beams  is  not  part  of  the  impulse  response,  we  must  assume 

that  their  lengths,  and  ,  are  not  necessary  equals.  This  pre¬ 

sents  a  first  obstacle  to  the  adaptation  of  Levinson's  (1949)  classical 
proof  of  the  uniqueness  of  the  inverse  Strum-Lionville  to  the  problem  at 

hand.  Indeed,  in  order  to  compare  the  two  beams  we  must  bring  them  to  some 

common  ground.  If  and  x^  are  coordinates  suitable  for  each  beam, 

then,  we  can  accomplish  our  aim  by  introducing  a  new  variable  Z  as  follows: 


2.8 


r 


and 


(2.17a) 


(2.17b) 


Note  tnat  as  and  vary  from  0  to  respectively  and 

,  c  varies  from  0  to  H  which  was  defined  in  (1.39).  This  is  a 
consequence  of  the  fact  that  eigenvalues  of  the  two  beams  have  the  same 
asymptotic  behavior. 

The  functions  5(x^)  and  c(x^)  defined  in  (2.17)  are  nondecreasing, 
and  hence  are  invertible: 


x 


m  ,  x(l) 


(e) 


x‘2>  =  x‘2)(5) 


(2.18) 


We  are  now  ready  to  define  the  following  strange  integral  of  a  hybrid  version 
of  the  Green's  function  a  given  in  (2.8),  namely 


+i» 


(2) 


Jti)  «  2W  f  ds  f  P(2)Cx(2)Cc))  Qtz*1)(x(2)(0,  x(1)(0;  -s)  h(;)  dc 


-l® 


...  (2.19) 


i  . 


■■■  ^  !  V 


2.9 


«*'- '  ••  >•  ■  •  >  UWv; 


■ 


where 
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Fjn.V)^  Ml?  •  [  ^ 

m=l  um  2/2  s 5  [p°HL(  ')P  o(1,(0) 


(sinh  /2  sa  =  -  sin  J1  s^  E } 


(2.23) 


where  1=1,2.  Consequently 


[r0)(L(Vf2 


[p(1)(  L(1))r 


r(1)(0)  ‘4 

PO)(0) 


rC2)(0)  1  4 


[r(2)(l(2))]h  [p(2)(L(-2))f2  p  C2i  (0) 


The  same  procedure  applied  to  P1J(0,  u2),  viz 


(2.24) 


[r(,)(L(1))]2  F,  IT  (1*^^ 

m=l  m  2s2  [p  ^  (l/1 ' )] 


(cosh2(  =■)  -  cos2(  =)>  * 


(2.25) 


yields 


[ro)(Lo))j*s  •  fcUVV  [ru)aU))f  ‘  l2'26> 


In  the  same  vein,  we  can  confront  the  product  representation  (1.19)  for 
W(0,  -s)  with  its  asymptotic  representation  (1.26),  viz. 


[r('V’>)]*  n  0^),,,  l°(1iw  y'!1w|V0(^>)]^ 

m-1  m  I  p  (L  ] 


•  {cosh2  E  +  wS2  S)  , 


(2.27) 


from  which  we  deduce  that 


*  r.(i|(oi 


Making  use  of  (2.26)  we  can  rewrite  (2.24)  and  (2.28)  as  follows 

r(1)(0)  ,.(2)(0) 

p'^ST  =  ’ 

rO)(0)  p(1)(0)  =  r(2)(0)  p(2)(0)  , 

which  imply  (2.21)  and  (2.22).  Therefore,  if  the  two  beams  clamped  at  the 
right  end  have  the  same  impulse  response,  their  densities  and  rigidities 
must  coincide  at  the  left  end.  Therefore,  we  are  justified  in  omitting  the 
superscripts. 

In  addition,  they  must  satisfy  (2.26),  i.e. 

rO)(L0))pU)(L(l)),r(2)(|_(2))o(2)a(2))  . 


Note  that  the  above  product  is  finite  even  for  the  pathological  beams  pre¬ 
viously  discussed. 


2.12 


Evaluation  of  J(z)  by  Calculus  of  Residues:  An  examination  of  the 
integrand  of  Jl  ?)  reveals  that  the  path  of  integration  can  be  closed  by 
a  half-circle  in  the  left  half  of  the  s-plane.  The  only  singularities 
are  simple  poles  at  s  =  and  a  straight-forward  application  of  the 
calculus  of  residues  yields: 


,(')(xll) 


“S' 


»<2>(*<2)> 


5(2>u(2))h 


dx 

dc 


(2) 


dt 


p  (2)^(2)}  ^2)(xt2))h 


dx^ 

dS 


dt 


> 


trC2)(L(2))]2n  (lA2-) 


m^n 


m  J 


...  (2.29) 

In  the  above  formula,  we  have  not  bothered  to  write  explicitly  that  x^ 

(2) 

and  xv  '  are  functions  of  5  and  c  .  We  have  also  used  the  representation 
(1.19)  for  the  Wronskians. 

At  this  stage,  we  make  use  of  the  fact  that  ^^(x^)  and 
^ (x^1 for  1=1,2  are  related  by  a  constant  of  proportionality  kn  , 
which,  as  (1.32)  reveals,  is  soley  dependent  on  the  impulse  response  data,  and 

hence  identical  for  both  beams.  As  a  result  (2.29)  becomes 


2.13 


3U)  • 


J, 


n  (1  - 
mj<n 


n  \ 

^iL 
m 


[rt,!(LC11)j'2  ' 


-  /2)(xt2),  w2)]h  ^  }  dc;'r^(x{1)) 

+  [/2V2))r2  ;f  Pt2)oc(2))fJ2,(x(2))h42)  d;~;  • 

[kn  4>(1)(x(1),  w2)  -  ^(1)(x(1),  w2)] 


The  above  expression  can  be  further  simplified  by  means  of  (.1.34'),  i.e. 
by  recalling  that 


kn*(,)cxCi).^).»(i)(x(i)j„2)  =  km*,(i)Cxm) 


where,  as  (1.36)  indicates 


[r0)(L0))]2  kO>*  ,  [r(2)(L(2))]2  k(2)* 


n 


Fi  n(l  -  ) 

m  m 


As  a  result 


2.14 


3(0  'Jl  CL^l? 


Wn 

n  n 


pC2)^xC2)j  /2)(xt2))h  ^  dc 

£o _ _ _ ^ _ 

2  2 
OJ  oj 

n  0  -  -7-  )  nCl  -  -r  ) 
m^n  “mm  in 


(2.30) 


To  avoid  problems  in  the  pathological  cases,  it  is  preferable  to  free  the 
above  expression  of  the  terms  r^(L^)  r^(L^)  .  This  can  be  done 
by  defining  the  norm  of  ,  viz. 


y(i)i  p 

n  1  1 


)[*n(i)(x(1) 


)]2 


(2.31) 


The  superscripts  are  not  necessary  for  the  actual  evaluation  of  the  above 
norms;  hence  they  shall  be  omitted  temporarily.  We  recall  the  equations  for 
4 (x ,  -s)  and  ?n(x),  viz. 


and 


d2  „  d24 
dx2  r  dx2" 


-S  p  4  , 


d2 

dx2 


d2  fr 
~dP~ 


-  w2  p  T 

n  M  n 


Multiplying  the  first  by  Tn  ,  the  second  by  4  ,  subtracting,  inte¬ 
grating  the  resulting  expression  over  (0,  L),  we  get  after  making  use  of 


2.15 


0.22)  and  (1.23) 


,(0)  r(0)  *,,,(0.  s)  =  u£(l  +  V  fL  p  *  i 

•* 


Tn  dx  • 


(2.32) 


The  product  representations  (1.28)  and  (1.30b),  enable  us  to  write: 


a  (1  +  T-)  n(1  "  Sr)  =  fL 


n  m^n  m  m 


p  *  v  dx 

n 


It  only  remains  now  to  let  s  •+  -w*  to  see  that 


ri^Ui  i in  -1)  id  -  $>  =  k„  ll«nll2 

n  m^n  m  m 


(2.33) 


and  consequently  (2.30)  becomes 


[  p^(x^)  f^U^h  ^  dc 

1 0 _  n  ^ 

ill’ll  H^’l! 


3(0  >  l 

n=l 


i1}(x0))  .  (2.34) 


Needless  to  say,  if  the  two  beams  were  identical,  then  we  could  drop  the 
superscripts  altogether  and  get  the  classical  eigenfunction  expression 


JtiW)  =  I 

n=l 


|  p  %  h(s(x))  dx 


fn(x)  , 


2.1 


as  a  special  case  of  that  more  general  eigenfunction  expansion.  Incidentally, 
in  this  particular  case,  JCc)  is  just  hU)  . 

Evaluation  of  3(0  via  residue  of  pole  at  infinity:  An  alternative 
expression  for  is  obtained  by  taking  the  contour  in  the  s-plane  to  be 

a  very  large  circle  and  replacing  ^(x^,  x^;  -s)  by  its  asymptotic 

expression  for  large  (s|  .  The  calculation  is  straightforward  by  tedious. 
Once  this  is  done,  we  perform  the  integration  over  c-from  (0,  5)  and 
U,  e)  asymptotically:  this  requires  a  mere  integration  by  parts,  provided, 
of  course,  that  r,  p  are  differentiable.  The  remaining  integral  over  s  is 
now  trivial  since  to  leading  order  the  s-dependence  of  the  integrand  is  1  . 
The  upshot  of  these  calculations  is 


Equating  the  two  expressions  for  f  (s)  we  get: 


r(2)(x(2)) 


f  „<2>  *<2>  h  §l2dt 
y  jo  ,  "  *0)(x0h 


...  C2.36) 


But  the  left  hand  side  can  also  be  written  in  terms  of  a  standard  eigen¬ 
function  expansion,  viz. 


2.17 


We  must  again  digress  before  we  can  draw  any  conclusion  from  (2.39). 


The  Liouville  transformation  and  the  canonical  4tn  order  equation:  Willy  nilly, 
we  have  been  pushed  first  into  the  introduction  of  the  variable 


(•A  _  . 

c  -  c  £  r*  dx 


and  now  into  considering  the  functions 


i  I  v  (x(0) 

unCO  =  Cptx(O)]  CrCxCO)] 


(2.40) 


We  can  think  of  these  new  variables  as  constituting  a  Liouville  trans¬ 
formation.  In  terms  of  them,  the  original  eigenvalue  problem,  viz. 


d*  rd2,r 
dx2  dx2 


w  2  p 
n  n 


■  (r  )'  ■  0  ,  at  x  =  0  , 


%  ■  ■ 0  • 


at  x  =  L  , 


is  transformed  into 


d4U  a  d  U 

T#  +  a?  Al5>  SC  *  u„  ■  “n  “n 


(2.41) 


+  b  o  =  o 


at  5  =  0  (2.42a) 


+  e  Un  =  0 


2.19 


and 


2.20 


Since  we  have  been  led  to  a  consideration  of  the  canonical  4^  order 
equation  (2.41),  one  may  ask  why  did  we  not  consider  this  equation  from  the 
beginning?  As  a  matter  of  fact,  I  have  considered  the  inverse  eigenvalue 
problem  for  this  equation  in  a  previous  paper  (Barcilon,  1974).  That 
analysis  was  valid  only  in  the  case  where  the  eigenvalues  were  simple. 

This  assumption  is  correct  provided  that  the  canonical  equation  can  be 
written  as 


d  d  d 
a°  d7ai  dTa2  d7as 


o4  U  =  U)  U 


(2.47) 


where  aQ,  a*,  a2  ,  a3  and  ai+  are  positive  and  if  the  boundary  conditions 
satisfy  certain  requirements  (Karlin,  1971).  Thus,  one  is  not  dealing  with 

1L 

the  most  general  canonical  4  order  equation,  but  with  that  subclass  which 
has  simple  eigenvalues.  Under  those  circumstances,  it  is  preferrable  to 
deal  with  the  beam  equation. 


Before  resuming  our  discussion  about  uniqueness,  we  shall  introduce 
some  new  variables  which  will  simplify  the  formulas  for  A(0,  B  (c ) ,  a,  ....  e, 
namely 


P(0 


& 

D  * 


(2.48) 


or  equivalently 


o(C)  =  o (0)  exp  ( 


P(s)  dc)  , 

0 


q(0  =  q(0)  exp( 


Q(c)  dc)  . 


(2.48') 


2.21 


With  these  variables,  (2.43)  becomes 


4Q5  +  -  2Q2  -  2QP  -  P2  -  A  =  0  (2.49a) 


+  Qp«  +  ‘  ^  + 

(2.49b) 

-  4QQ?P  +  Q?P2  +  2QPP?  -  2Q2P, 

+  Q4  +  2Q3P  -  4Q2P2  -  B  =  0 


or  better  still 


where 


Qr  =  F(Q)  +  A  (2.50) 


(2.51) 


F  = 


-  4Q?5  +  4QQS  +  2Q?P  +  2QP ^  +  2PP? 
+  4QQ5?  -  q5|  -  Q2  +  2QQ?P  -  4QPP5 
-  Q^P  +  Q4  +  2Q3P  -  4Q2P2 


(2.52) 


2.22 


i 


Finally,  the  equations  (2.46)  can  be  written  as  follows 


Q  (0 )  = 

% 


c  "I 

a  -  2c 

b  +  c2  -  ac  | 

d  -  3b  -  2c2  +  ac 
e  +  cb  -  ac2  -  cd  +  2c3 


(2.54 


Uniqueness  results:  Let  us  summarize  the  current  state  of  affairs.  Given 
two  beams  p^(x^),  r^(x^)  and  p^(x^),  r^(x^ )  with  the  same 
impulse  response,  then,  provided  that  these  functions  are  sufficiently 
smooth,  we  can  make  a  Liouville  transformation  and  write 


The  fact  that  Un  has  no  superscript  is  a  consequence  of  (2.39),  viz. 


qU)(0  ||fnC1)||  UnC5)  qC2](5)  jjy(2) 


for  every  n. 


The  coefficients  A^.B^  and  A^,B^  are  obtained  by  substituting 
p^,  q^  and  p^,  q^  in  (2.43).  Now,  if  we  were  to  subtract  (2.55b) 


from  (2.55a)  we  would  get 


(A^  -  A^)  (B^  -  B^ )  Un  =  0  . 


(2.56) 


Since  Un  satisfies  (2.42b),  we  are  forced  to  conclude  that 

A^U)  =  A^U)  *  A(?)  , 

B(1,(0  =  B(2)(0  -  BU)  • 


(2.57) 


or  else,  Un  would  be  identically  zero.  In  a  similar  way,  a  consideration 
of  the  boundary  conditions  (2.42a)  at  5=0,  implies  that 


a<»  •  .(«  ■  ,  , 
b<’>  =  b<2>  =  b  , 

C<-1)  =  =  c  . 

d<’l .  d«>  =  d . 
e<»  ,  «(2)  ,  e  . 


(2.58) 


(This  can  be  seen,  for  instance,  by  considering  eigenfunctions  associated 
with  large  values  of  n  .)  As  a  result,  the  vectors  and  are 
solutions  of  the  same  differential  equation  (2.50)  with  initial  conditions 


2.24 


(2.54).  Assuming  that  the  solution  of  this  initial  value  problem  exits, 
it  is  unique  since  F  is  Lipchitz.  Consequently 


and  in  particular 


Q(1)(£)  «  Q(2)(£) 

%  ^ 


Q(1)(£)  -  o(2)U)  , 
pO)(5)  =  p(2)(c)  . 


and  finally,  in  view  of  (2.21)  -  (2.22) 


o(1)(x(1)U))  =  o(2)(x(2)(0)  , 
r(1)(x(1)(c))  =  r(2)(x(2)(0)  • 


(2.59) 


(2.60) 


For  the  smooth  class  of  functions  which  we  have  been  considering,  we  can 
in  fact  go  further  and  see  that  ,  and  hence 


2.25 


1*1 


0(1)(x)  =  P(2)(x)  ,  1 

r0)(x)  =  r(2)(x}  . 


(2.61) 


However,  because  of  the  pathological  beams  previously  discussed,  I  suspect 

f 

that  (2.60)  is  the  best  result  possible.  Thus,  if  we  are  willing  to  dis-  /j 

regard  the  pathological  beams,  we  can  say  that  the  impulse  response  deter¬ 
mines  a  beam  uniquely. 


EXISTENCE  AND  CONSTRUCTION  OF  SOLUTION 


In  this  section  we  shall  outline  two  procedures  for  reconstructing  a 

beam  from  its  impulse  response.  One  of  these  procedures  was  previously 

presented  by  Barcilon  (1979a,  1979b).  The  central  idea  revolves  around  the 

use  of  continued  fractions  and  owes  much  to  the  work  of  Krein  (1951,  1952a). 

The  procedure  for  reconstructing  c (x)  and  r(x) ,  given  F  ,  F?  and 

(w  ,v  ,y  }  ,  will  work  provided  that  the  data  are  bona  v Lde.  data.  Thus, 

n  n  n  1 

we  shall  have  to  elaborate  criteria  for  recognizing  whether  or  not  a  solution 
exists  for  a  given  impulse  response.  We  do  not  have  definitive  results  about 
all  such  criteria:  we  shall  derive  a  few  and  suggest  the  need  for  more. 

Before  embarking  on  this  program,  it  might  be  useful  to  review  the 
situation  for  the  inverse  problem  for  a  vibrating  string.  The  impulse 
response  for  a  string  in  the  free/fixed  configuration  consists  of  (i)  the 

00  DO 

length  L  of  the  string  and  (ii)  the  spectra  { X  } ^  and  {y  associated  with 

the  fixed/fixed  and  free/fixed  configurations.  The  numbers  L,  {y  ,\ 

n  n  1 

constitute  a  bona,  ^idz  impulse  response  if  they  satisfy  the  following 
conditions : 

(i)  asymptotic  behavior 

a  HIT  , 

A  *\*  y  ^  — - -  ,  as  a  ®  (j. 

n  n  Li  ^ 

/  p  2(x)dx 

0 


(ii)  interlacing 


(iii)  global  condition 


oo 


The  asymptotic  condition  is  easy  to  understand  and  can  be  deduced  in  a 
straightforward  manner  form  the  direct  problem.  The  interlacing  condition 
is  also  sensible:  after  all  the  u's  are  associated  with  a  configuration 
freer  than  that  for  the  \'s.  Less  obvious  is  the  fact  that  this  condition 
insures  the  positivity  of  p (x) .  This  is  connected  to  a  theorem  of 
Stieltjes  (1894)  which  we  shall  discuss  in  the  sequel.  The  third 
condition  insures  that  the  total  mass  of  the  string  is  finite.  It 
plays  an  essential  role  in  the  construction  of  the  solution  to  the 
inverse  problem  for  a  vibrating  string.  We  would  like  to  find  the 
conditions  for  the  beam  analogous  to  (3.1)  -  (3.3)  for  the  string. 


Interlacing  of  eigenvalues. 


One  would  expect  that  the  following  ordering 


would  hold: 


w,  <o.  <v,  <u,  <  X,  <  ...  <  U).  <o.  <v.  <u.  <  X.  <  ... 
11111  11111 


However,  this  is  not  the  case.  Indeed,  by  considering  the  homogeneous 
beam  (i.e.  o  and  r  constants),  we  can  see  that  the  spectrum  associated 

with  the  clamped/clamped  configuration,  namely  ,  does  not  inter- 
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lace  with  that  for  the  f ree/claroped  one,  i.e.  {u)  ^  (Platzman, 
private  communication) . 


3.2 


To  elucidate  this  question,  we  appeal  to  a  paper  of  Krein  (1939). 
If  we  were  to  apply  the  theorem  stated  in  that  paper  to  the  case  of  the 
beam  at  hand,  we  would  get: 


w.  <  a.  <  .jj  ,  (3.4a) 

1  i  i+2 


a  .  <  X  .  < 

i  l 


i+2  ’ 


(3.4b) 


v  <  A  < 

l  l 


i+2  ’ 


(3.4c) 


u  <  A  <  u 

l  i  i+2 


(3 . 4d) 


These  results  are  general  results  stemming  from  the  properties  of  Green's 
functions  of  fourth  order  operators  which  are  so-called  "oscillating". 

These  results  can  be  improved  for  the  particular  case  of  the  beam  operator. 
The  first  step  consists  in  examining  the  first  eigenvalues  of  our  five 
eigenvalue  problems.  We  start  by  considering  and  To  that  effect 

we  introduce  a  hybrid  eigenvalue  problem  as  follows: 


(rv")"  =A^py,  0<x<L  , 

ay-(l-ct)y'  =  y"  =  0,  at  x  =  0  , 

y  —  y  *  *  0 ,  at  x  =  L 

where  a  is  a  parameter  ranging  over  (0,1). 

Clearly 


(3.5) 


(3 .  t>a) 


3.3 


and 


After  differentiation  with  respect  to  the  parameter  a,  we  can  easily 
see  that 


a(l-a)  /  pvndx 
0 


We  shall  specialize  this  result  to  the  first  eigensolutions  of  (3,5) 
and  (3.7)  and  appeal  to  a  result  of  Gantmakher*  (1936),  namely  that 


y1(x,a)  >  0  , 

for  0  <  x  <  L 

n1(x,a)  >  0  , 


(3.9) 


*Strictly  speaking,  we  should  prove  first  that  (3.5)  and  (3.7) 
can  be  transformed  into  integral  equations  with  oscillating  kernels. 
This  is  indeed  the  case  thanks  to  a  theorem  of  Karlin  (1971). 


3.4 


Consequently 


!  . 


yx  (0,a)  _>  0 


and  in  view  of  the  boundary  condition  at  x=0  , 


yi’  (0,a)  i  0  . 


We  prove  next  that 

(rn,  )  (0,a)  _<  0  . 


(3.10) 


(3.11) 


Let  us  assume  the  contrary.  Then,  by  continuity  there  exists  an  interval, 
say,  (0,5),  over  which 

II 

rnL  >  0  . 

This  interval  must  be  smaller  than  the  length  of  the  beam.  Indeed,  other¬ 
wise  would  be  an  increasing  function  throughout  (0,L)  and  since 
n^’(0,a)  must  be  positive,  It  would  follow  that  n^’(L,a)  is  also  positive. 
But  this  is  impossible  on  account  of  the  boundary  condition.  Thus 
0  <  5  <  L.  We  can  be  more  specific  and  define  5  such  that 

(rn^  )(5,a)  *  0  . 


Now,  integrating  (3.7)  twice  from  0  to  5,  we  see  that 


and  since  (rn^  )(0,a)  is  positive,  it  follows  from  the  boundary  condition 

»r  * 

that  (rp^  )  (0,a)  is  also  positive.  Hence,  the  left  hand  side  of  the  above 

equation  is  negative  whereas  the  right  hand  side  is  positive.  Therefore 
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(3.11)  must  be  true  and  as  a  result,  the  derivative  of  with  respect 
to  a  is  positive,  i.e. 


V1  "  U1 


(3.12) 


By  considering  a  problem  similar  to  (3.5)  but  with  the  boundary 


conditions 


ay  -  (l-a)y’  =  (ry")  =0  at  x  =  0 


we  could  show  that 


°1  "  V1 


(3.13) 


In  a  similar  vein,  we  could  also  prove  that 


W1  <  °1 


(3.14) 


Combining  all  of  these  results,  we  get 


<  °1  <  V1  <  U1  <  X1 


(3.15) 


The  technique  used  to  establish  (3.15)  is  not  suitable  for  the 
higher  eigenvalues.  To  order  them  I  have  found  that  the  introduction 


of  certain  auxiliary  functions  is  very  helpful. 


3.6 


Auxiliary  functions.  Recall  that  u(x,-s)  and  v(x,-s)  were  fundamental 
solutions  of  the  beam  equation  chosen  in  such  a  way  as  to  satisfy  the 
clamped  boundary  conditions  at  x=L  (see  (1.7)).  Recall  also  the  forced 
problem  (2 .3)— (2- 4)  which  led  to  the  impulse  response  when  p(x)n(x)  was 
set  equal  to  6(x).  We  can  rewrite  this  problem  as  follows: 


II 

(ry")  *  -spy  , 


ry' 


(ry”)  -  1 


at  x=*0 


(3.16) 


y  *  y’  *  0,  at  x=L 


We  can  avail  ourselves  of  the  functions  u  and  v  to  solve  this  problem. 
In  fact 


y(x,-s)= 


u(x,-s) (rv”) (O.-s)-v(x.-s) (ru”) (0,-s) 

(ru”) (0,-s) (rv")  (G,-s)-(rv') (0,-s) (ru")  (0,-s) 


and  as  a  result,  the  impulse  response  given  in  (2.9)  can  also  be 
written  as 


(3.17) 


u(0,-s) (rv") (0,-s)-v(0,-s) (ru") (0,-s) 

y(0,-s)=  —  ■■■  ■  - — - 

(ru") (0,-s) (rv") ' (0, -s)-(rv") (0,-s) (ru") ’ (0,-s) 


(3.18a) 


and 


9(0, -s)« 


u’ (0,-s) (rv") (0,-s)-v’ (0,-s) (ru") (0,  -s) 
(ru") (0,-s) (rv") ' (O.-s)-(rv") (0,-s) (ru") ’ (0,-s) 


(3.18b) 


3.  7 


These  expressions  suggest  that  we  define  the  following  functions: 


Y(x,to2}  =  u(x,u>2)  (rv")  (x,u>2)  -  v(x,w2)( ru”)(x,^2)  , 


(3.19) 


e(x,J)  =  u’  (x,uj2)  (rv")  (x,uj”)  -  v’ (x,o)2)  (ru")  (x,ai2)  ,  (3.20) 


D(x,a)2)  =  (ru”)  (x,u)2)  (rv")  ’  (x,u>2)  -  (rv")(x,.2)(ru”)’(x^2),  (3.21) 


and  that  we  work  directly  with  them  rather  than  with  u  and  v.  Numerically, 
this  approach  has  the  advantage  of  minimizing  the  loss  of  accuracy  due  to 
cancellations  in  the  computation  of  Y,  0  and  D. 

We  can  look  upon  Y,  0  and  D  as  determinants  obtained  by  con¬ 
sidering  the  2x2  submatrices  arising  from  the  4x2  matrix 


ru  rv 


(ru")’  (rv")’1 


Three  other  such  determinants  can  be  formed,  namely 


I(x,u)2) 


(3.22) 


J(x,00  ) 


u  v 


(ru")’  (rv")’ 


(3.23) 


3.3 


and 


K(x,uj  ) 


(ru")’  (rv")' 


(3.24) 


It  is  a  simple  matter  to  check  that  these  six  auxiliary  functions 
satisfy  the  following  differential  equations: 


I-iv  . 

(3.25a) 

Y’  =  0  +  J  , 

(3.25b) 

J’  *  K  , 

(3.25c) 

0’  -  K  , 

(3 . 25d) 

K’  =  i  D  -  W2pl  , 

(3.25e) 

D’  =  -io2dY  , 

(3. 25f ) 

where  primes  indicate  differentiation  with  respect  to  x.  These 
auxiliary  equations,  rather  than  the  beam  equation,  will  occupy 
the  center  stage  during  the  reconstruction  procedure.  Note  that  in 
view  of  (1.7),  the  boundary  conditions  associated  with  (3.25)  are 

I*Y*J  =  0  =  K  =  D-rZ(L)  =0  at  x-L.  (3.26) 


Keeping  in  mind  the  definitions  (3.19)— (3.24)  of  the  auxiliary 
functions,  it  is  a  simple  matter  to  show  that 

00  ^ 

KO, J)  -  r2(L)(FoF2-F12)  ]  ["  (1-  iiV-  )  , 

n=l  ^n 


(3.27a) 


3.9 


In  the  above  expressions  and  F2  are  the  first  and  second  moments  of 
the  flaccidity  defined  in  (1.27)  and  (1.29)  while  Fq  is  the  zeroth 
moment ,  namely 


L 


0 


(3.28) 


Note  that  (3.25c)  and  (3.25d)  together  with  the  boundary  conditions 
(3.26)  imply  that 

0(x,u)2)  =  J(x,uj2)  .  (3.29) 


This  result  is  easy  to  understand:  it  means  that  the  eigenf requencies 
of  the  truncated  beam  (x,L)  are  identical  if  the  boundary  conditions 
at  the  left  end  are  y  »  (ry") '  *  0  or  y’  =  ry"  =0.  As  we  already 


3. 10 


know,  these  two  problems  are  adjoint  of  one  another. 
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In  deriving  the  expression  (3.27a)  for  1(0, ui  )  we  made  use  of  an 
identity  for  the  auxiliary  functions  which  will  play  a  crucial  role  in 
the  sequel,  namely 


-I (x,u)2)D(x,uj2)  +  Y(x,u)2)K(x,uj2)  _  02(x,w2)  =  0 


(3.30) 


This  identity  can  be  established  from  the  differential  equations  (3.25). 
However,  its  character  is  algebraic  in  nature  and  it  is  best  obtained 
from  the  very  definitions  of  the  auxiliary  functions. 

The  product  representations  (3.27)  for  x=0,  can  be  viewed  as 

special  cases  of  product  representations  for  arbitrary  x's.  Indeed, 

2 

the  auxiliary  functions  are  entire  functions  of  w  of  order  1/4.  Con¬ 
sequently,  we  can  write 


I(x,u2)«r2(L)[  | 


„  2,  2 

“dt  1  (1 - f -  ), 

r  J  1  A  “(x) 


(3.27 ’a) 


Y(x,w  )  -  -r  (L) 


rL  2  ■ 
]  ^  « 


u  2(x) 
n 


(3.27’b) 


0 (x,ui2)  -  r2 (L) 


v  "(x) 
n 


(3.27'c) 


3.11 


3.12 


T 


Indicating  a  derivative  with  respect  to  a  by  a  dot,  we  deduce  that 


(.ry  ")"  =  2X  X  py  +  X  py 
n  n  n  Jn  n  n 


(3.31) 


and  from  the  boundary  conditions: 


y  +y’=y’+y  =0  at  x  =  a  , 

7n  Jn  7n  7n 


y  =y’=0  at  x  =  L  . 
'n  n 


Multiplying  (3.31)  by  y^  and  integrating  over  (a,L)  we  get: 


-r(a)yn"(a,a)yn’  (a, a)  =  2XnxJlynii 


which,  after  some  simple  manipulations,  yields 


dXfl(a)  r(a)[yn  (a, a)]' 
^  2Xn(a)  2 


A  similar  procedure  could  be  used  for  the  other  spectra. 


The  grand  Interlacing.  We  are  now  in  a  position  to  resume  our  discussion 


about  the  interlacing  of  the  various  spectra. 

2  2 

Let  us  consider  the  curve  I(x,w  )  for  a  fixed  value  of  w  lying 


between  X 


^  and  X £2  (see  Fig.  2). 


Insert  Fig. 2 
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I(x,u>  )  is  positive  for  x  in  (0,1)  and  negative  for  x  in  (2,L)  where 


2  ? 
a)  =  </. 


2  2 

(3.25a)  implies  that  the  zeros  of  I’  coincide  with  the  curves  u  =un_(x). 
Therefore  I*  vanishes  only  once  in  (8,,L)  since  an  immediate  generalization 
of  (3.15)  would  yield 


u  (x)  <  A1(x)  . 


Over  (0,8.),  I*  could  either  vanish  or  not  vanish.  If  I’  were  different 

2 

from  zero  for  all  x's  in  (0,8.(w  ))  and  for  all  id's  in  the  interval  (A^,A0), 
then  we  would  have 


u2  >  *2 


which  contradicts  (3.4d).  On  the  other  hand  if  I*  vanished  three  or 
more  times,  then  we  would  have 


y4  <  X2 


which  also  contradicts  (3.4d).  Therefore  I*  can  vanish  either  once  or 
twice  in  (0,8.).  Let  us  examine  the  second  case  more  closely.  Graphically, 
the  situation  is  as  indicated  in  Fig.  3  .  Consider  the 


Insert  Fig.  3 


the  structure  of  I  on  the  straight  line  u>  =  A^.  Certainly 
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(3.32) 


1(0, x2  )  =  0  , 


and 


I  (x ,  A  )  <  0 


(3.33) 


over  the  interval  (0,A(Ao  )).  Now 


Y(0,A2  )  >  0 


But,  because  of  the  auxiliary  equation  (3.25a),  this  inequality  implies 
that 


I’ (0,A2  )  >  0 


and  hence  for  small  values  of  x,  it  follows  from  (3.32)  that 


I(x,A7^)  >  0 


which  contradicts  (3.33).  Hence  I*  can  only  vanish  once  and  as  a  result 


W1  < 


A ,  <  u  < 
1  u2 


Repeating  this  procedure,  we  can  conclude  that 


Ui  < 


Xi  *  Ui+1  " 


i+1 


(3.34) 


3.15 


- .  -  —  *  •  V«t '.a. •  **r  ■  - -»  . ►-ri|y. 


By  exploiting  (3.25b) 


or  rather 


Y’  =  20  , 


we  could  show  in  a  similar  manner  that 


j .  < 

1 


ui 


v  . 


+1 


‘  i+1 


Similarly,  (3.25d)  and  (3.25f)  would  yield 


a .  < 

l 


v . 

x 


<  a 


i+1 


i+1 


(3.35) 


(3.36) 


and 


“i  "  yi  " 


i+1 


< 


i+1 


(3.37) 


Additional  interlacings  can  be  obtained  by  means  of  the  quadratic 
identity  (3.30).  Indeed,  setting  x=0  and  u=wn  in  (3.30),  we  see  that 

Y (0,w  2)  K(0,u  2)  -  02(0,o)  2)  =  0  , 
n  n  n 

which,  together  with  (3.37)  implies  that 


“i  " 


°i  " 


“i+1 


"  °i+l 


Had  we  set  x=0  and  in  (3.30),  we  would  have  obtained 

-I(0,o  2)  D(0,o  2)  -  02(O,o  2)  -  0  , 
n  n  n 


(3.38) 
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which,  together  with  (3.38)  implies  that 


o,  <  X .  < 
i  1 


i+1 


i+1 


(3.39) 


By  combining  all  of  these  inter lacings  we  can  write  the  srand  inter¬ 
lacing 


<  a .  < 

i 


vi  "  ui  " 


X, 

( -U) 


°i+l 


V  <  u  < 

i+1  Mi+1 


l  i+l') 

K+J 


(3.40) 


The  bracket  (  u)*+^  )  indicates  that  the  order  between  these  two  eigen¬ 
values  cannot  be  decided.  As  we  alluled  to,  the  homogeneous  beam 
provides  an  instance  where  one  can  check  that  the  X’s  and  the  u>’s  do  not 
interlace. 

Clearly,  the  data  ought  to  satisfy  (3.40)  in  order  for  a  solution 

to  exist.  More  specifically,  since  only  { cu  },  {v  }  and  }  are  given, 

-  n  n  n 

the  data  should  satisfy  the  small  interlacing,  namely 


“l  <  V1  <  U1  <  “2  <  V2  <  W2  < 


(3.41) 


We  are  prepared,  of  course,  to  require  that  the  ui,  v  and  u  spectra  have 
the  asymptotic  behavior  given  in  (1.38)  and  perhaps  some  other  gross 
condition(s)  akin  to  (3.3)  for  the  vibrating  string.  But,  as  we  shall 
soon  see,  the  situation  is  far  more  complicated  and  other  conditions 
must  come  into  play.  Some  of  these  additional  conditions  arise  quite 
naturally  in  the  process  of  constructing  the  solution  to  the  inverse 
problem.  Therefore,  we  must  examine  this  process;  the  first  step  is  a 
discretization  of  the  direct  problem. 
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A  useful  discretization.  The  basic  idea  for  solving  the  inverse  problem. 


is  to  look,  upon  the  beam  with  given  spectra  {ou  ,v  ,  p  }.  as  the  limit 

n  n  n  1 

N-l 

of  a  simpler  beam  with  truncated  spectra  {w  ,  v  ,  u  ),  as  the 

n  n  n  I 

number  of  eigenf requencies  tends  to  infinity. 

Beams  with  a  finite  number  of  eigenfrequencies  must  have  a 
finite  number  of  degrees  of  freedom.  Thus,  they  must  be  made  up  by  a 
finite  number  of  point  masses.  Consequently,  their  density  structure  is 
of  the  form 


N-l 

p (x)  =  £  m.  6(x-x.). 

i=l  1  1 

The  structure  of  the  flexural  rigidity  is  not  dictated  as  clearly.  One 
useful  approximation  is 


1 

r  (x) 


i=l 


f .  5 (x-x . ) 
i  l 


I  have  discussed  already  such  discrete  N-  beams  (Barcilon  1979a, b) : 
they  can  be  thought  of  as  being  made  of  N+l  weightless  segments  of 

length  { £ ^ } o  and  of  infinite  rigidity,  connected  by  clothespin-like 

N— 1  N 

devices  of  mass  {m^}^  and  "flaccidity"  (or  limpness)  (f^l^  (see 

Fig.  A  ). 

Insert  Fig.  A 


These  elastic  joints  are  located  at 


i  ,N 

{Vi 


where 
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(3.42) 


(3.43) 


xi+l  -  xi  =  Z1  ’  (i“0’  ••••  N) 

It  is  important  to  realize  that  the  structure  (3.42)  of  the  density  is 
solely  responsible  for  the  finite  number  of  eigenf requencies  of  the 
resulting  beam.  The  structure  (3.43)  of  the  flaccidity  is  chosen  for 
convenience.  In  fact,  other  forms  could  have  been  chosen  which  might 
indeed  be  preferable  in  other  problems.  For  instance,  we  could  have 
written 


1 

r  (x) 


for 


x.  ,  <  x  < 
l-l 


X  . 
1 


i.e.  we  could  have  considered  the  weightless  segments  to  have  finite 
constant  rigidity.  This  case  is  briefly  discussed  in  an  appendix. 

Since  we  are  interested  in  the  discrete  version  of  (3.25)  it 
would  seem  expedient  to  substitute  in  these  equations  p (x)  and  r(x)  by 
their  expressions  as  given  in  (3.42)  and  (3.43).  However,  this 
approach  is  dangerous  because  of  the  occurrence  of  products  of 
generalized  functions.  We  shall  follow  therefore  the  long  route 
and  discretize  (1.2)  first. 

Over  a  generic  interval  (x^^,  x^)  the  equations  (1.2)  become 


y*  =  Q  . 

9’  =  0  , 
t’  =  -X> 

X*  =  0  . 

Integrating  once,  we  see  that  for  x^  <  x^  , 
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(3.44) 


(3.45) 


t  'MS 


y(x,w2)  =  ei_1(w2)  [x-x.^]  +  Yi_1(“2) 


0  (x,0>2)  =  0i_1  (U>2)  , 


2  2  2 
t(x,io  )  =  )[x-x  ]  +  1^(0)  ) 


2  2 
X(x,4j  )  =  X1_1(w  )  • 


Clearly  y^  9^  Ti_^  anc*  X^  ^  correspond  to  the  values  of  v,  9, 

t  and  x  at  x^  ^+0.  At  x=x^,  the  following  jump  conditions  hold: 


yCx^+O.w2)  -  y (x1-0,ui2)  =  0  , 


9(xt+0,u2)  -  SCx^O.o)2)  =  f ±t ±  , 


2  2 

rCx^+OjiD  )  -  t(Xj-0,w  )  =  0  , 


2  2  2 
X  (x^+0, u)  )  -  xCx^^-O.u  )  =  -m^u)  y 


Substituting  the  expressions  (3.46)  in  the  above  jump  conditions  we 
obtain  the  discrete  version  of  (1.2),  namely 


yi  "  yi-l  +  li-l  9i-l  » 


9i  "  9i-l  +  fi  Ti  ’ 


Ti  “  Ti-1  "  *i-l  xi-l  ’ 


x i  *  *i-l  ~  -  mi  yi 


(3.46) 


(3.47) 


(3.48) 
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k 


We  next  define  two  fundamental  solutions  of  the  above  difference 


2  2 

equations,  closely  related  to  the  functions  u(x,id  )  and  v(x,u>  ). 

We  denote  these  solutions  by  a  superscript  (3)  and  (4),  and  characterize 
them  by  their  values  at  x=x^+^: 


= 

9(3)  = 

t(3) 

N+l 

N+l 

N+l 

= 

9^)  = 

t(4) 

N+l 

N+l 

tN+1 

i  -  <3)  _  n 
1  XN+l  ~  0  ’ 


+  1  =*  0  , 


or,  in  view  of  (3.46) 


(3)  =  (3)  „  (3)  _  =  (3)  = 

yN  N  N  1  XN  U  ’ 


y^4)  =  9<4> 
yN  °N 


+  *N  =  XN+1  =  ° 


Pursuing  the  analogy  with  the  continuous  case,  we  define 


I^u.2) 


43> 


,(3) 


yi 


(4) 

i 

(4) 


Y1  (uj2) 


.(3) 


(4) 


(4) 


J^u2) 


i  y 


.(3) 

i 

(3) 


.<*>  i 


(4) 


(3.49) 


(3.50) 


(3.51a) 


(3.51b) 


(3.51c) 
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e^u.2) 


,(3)  .(4) 


(3) 


0 


1 

.(4) 


K.  (.J) 


Di(J) 


(3) 

i 

e(4) 

i 

(3) 

,  (4) 
Xi 

i 

(3) 

x<4) 

i 

Ti 

(3) 

i 

(4) 

Xi 

(3 . 51d) 


(3 . 31e) 


(3 . 51f ) 


Making  use  of  (3.48),  we  can  deduce  that  these  discrete  auxiliary 
variables  satisfy  the  following  difference  equations: 


Xi  "  *1-1  +  fiYi  * 


(3.52a) 


Y  «  Y.  .  +  2..  .  0.  ,  +  i.  .  J.  , 

i  l-l  l-l  l-l  l-l  i 


Ji  =  Ji-1  +  4i-l  Ki-1  ’ 


(3.52c) 


®i  ®i-i  +  ^i-1  ^i-1  * 


(3 . 52d) 


Ki  "  Ki-1  +  fiDi  *  “  mi  ri-l  ’ 


(3. 52e) 


Di  =  Di-1  '  “  mi  Yi 


( 3 . 52  f ) 
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The  advantage  of  following  this  procedure  is  now  clear:  we  do  not  have  to 
worry  about  whether  the  derivatives  are  to  be  discretized  as  forward,  back¬ 
ward  or  centered  differences.  The  choice  is  made  quite  naturally.  As 
an  additional  bonus,  we  can  check  that  the  discrete  version  of  the 
quadratic  identity  (3.30),  namely 

-IiDi  +  YiKi  -  Q±2  -  0  (3.53) 

is  compatible  with  the  equations  (3.52).  Equation  (3.53)  takes  into  account 
the  face  that 


-  Qj  ,  (3.54) 

which  follows  from  (3.52c),  (3.52d)  and  the  end  conditions 

XN  "  YN  "  JN  "  ®N  "  S  =  V1  =  °  ‘  (3'55) 

Eventhough  the  generic  N-beam  under  discussion  has  a  length 
Xjj+^,  we  ended  up  applying  the  boundary  conditions  at  As  a  matter 

of  fact,  the  length  of  the  last  infinitely  rigid,  weightless  segment 
ijj  as  well  as  the  mass  of  the  clast  clothespin  m^  do  not  enter  into 
the  equations  of  the  discrete  beam.  This  is  a  consequence  of  the 
clamped  boundary  conditions  at  the  right  end.  Indeed,  the  last 
segment  remains  stationary  during  the  vibration.  Similarly,  the  last 
clothespin  flexes  but  does  not  move  up  or  down.  This  situation  is 
reminiscent  of  that  for  the  pathological  beams  previously  discussed. 
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at  x  « 


1 


Starting  from  the  values  of  the  discrete  auxiliary  variables 

2 

and  making  use  of  (3.52),  we  can  show  that  I^(w  )  is  a 

2  2  2  2 
polynomial  in  id  of  degree  N-i-2  whereas  0^(w  ),  K^(id  )  and 

2 

D^id  )  are  polynomials  of  degree  N-i-1.  In  particular,  at  the  point 
x*0  which  is  associated  with  i=0,  we  have 


N  N 


N-2 


V“2)  ] 


n-1 


(i-  , 

X 

n 


N-1 


V“2)  '  'ti 


(1-  s-j)  , 


n=l 


N-1 


So<“2>  -[  I  *lfil 


U-  , 


n*l 

N-1 


K  (w2) 

o 


(1-  “)  , 


n=l 


N-1 


D  (u2)  - 
o  l 

n»l 


(i-  — )  . 


These  expressions  are  the  discrete  analogues  of  (3.27).  We  are  now  in 

a  position  to  outline  a  program  for  solving  the  inverse  problem.  Given 

the  impulse  response,  i.e.  given  {oj^,  v^,  ant*  F2’  we 

attempt  to  construct  that  N-beam  whose  eigenf requencies  are 
N-1 

{id  ,  v  ,  u  },  and  such  that 
n  n  n  1 
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(3.56a) 


(3.56b) 


(3.56c) 


(3 . 56d) 


(3. 56e) 


i 


IAmUmis 


— 


N 


N 

I 


1 


x.2f  .= 
1  1 


(3.57) 


In  other  words,  given  the  3N-I  bits  of  data  (to  ,  v  ,  y  ^  ,  F, ,  F„ 

n  n  n  1  1  2. 

we  set  out  to  find  the  3N-1  unknowns  {x.}^,  {m. {f.}?  . 

l  1  i  1  l  1 

2  2 

Construction  of  I  (u  )  and  K  (w  )  .  With  the  given  data,  we  can  construct 
_ o _ o _ 

2  2  2 

the  polynomials  Yq(w  ),  0Q(<y  )  and  ).  By  means  of  the  polynomials 

and  the  quadratic  identity  (3.53)  for  i=0,  we  can  find  the  value  of  Iq 
at  N-l  points.  Indeed,  by  setting  w=Un  in  (3.53)  we  deduce  that 


I  (u  2) 

o  n 


0  2  (u  2) 

o  n 

D  (y  2) 
o  n 


,  n-l,  2 . N-l 


(3.58) 


2 

I  (<Ot  which  as  indicated  by  (3.56a)  is  a  polynomial  of  degree  N-2, 
is  therefore  completely  determined: 


N-l 

i„<“2)  -  -  l 

n*l 


g>n  > 

D  (»  2) 
o  n 


N-l 


Mn 

N-l 


k^n 


2 


2 


(3.59) 


Note  that  the  zeros  of  the  above  polynomial  are  real  and  positive.  Indeed, 

2  , N-l 

since  u>.  and  y.  interlace,  the  terms  in  the  sequence  (D  (y  );, 
li  o  n  1 
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•> 

alternate  In  sign.  As  a  result  the  sign  of  Iq(u  “)  alternates  with  n. 

2  (N)  N— 2 

We  shall  denote  the  zeros  of  I  (m  )  bv  {A  }  ,  the  superscript  N 

o  '  n  1 

being  used  to  remind  us  that  these  eigenvalues  are  associated  with 

the  N-beam.  We  shall  drop  this  superscript  whenever  there  is  no 

possible  source  of  confusion. 

We  have  just  outlined  a  procedure  for  finding  the  eigenf requencies 

of  the  N-beam  in  the  clamped/clamped  configuration.  The  fact  that  we  can 

find  these  eigenf requencies  is  not  surprising  since  F^,  and 
N-l 

{ u>  ,  v  ,  u  },  define  the  N-beam  completely.  What  is  surprising  is 
n  n  n  i 

that  we  can  do  it  without  having  to  first  solve  for  the  structure  of 

2 

the  N-beam!  In  the  same  way,  we  can  deduce  K  (u>  )  and  its  zeros.  We 

should  point  out  that  the  fact  that  we  were  successful  in  determining 
2  2 

Iq(u)  )  and  Kq (oj  )  is  related  to  the  way  in  which  the  given  spectra 
entered  into  the  quadratic  identity.  As  discussed  in  Barcilon  (1979b) 
three  spectra  and  two  gross  constants  a  sufficient  to  determine  a  non- 
pathologic  beam,  provided  that  these  spectra  are  sympathetic,  i.e. 
provided  that  they  are  such  as  to  yield  two  other  spectra  from  the 
quadratic  identity. 

N-l  (N)  N— 2 

Having  established  the  interlacing  of  {un}^  and  ).  , 

i.e. 


u 


< 

1 


< 


1(N)  < 
AN-2 


N-l 


we  can  establish  in  a  similar  manner  that 


< 


(3.62) 


.  (N) 

w-2 


(N) 

N-l 


To  these  interlacing  relations,  we  can  add  those  that  the  data  must 
satisfy,  i.e.  the  small  interlacing  (3.41),  or  rather 


"1  < 


U1  < 


N-l 


N-l 


N-l 


(3.63) 


However,  the  grand  interlacing  (3.40)  is  equivalent  to  nine  interlacing 
conditions;  we  are  missing  the  following  three: 


.00 


V1 " 


(N) 

<  o„  ;  <  v„  . 

I  iN— ± 


(3.64a) 


o 


(N) 

1 


“l  < 


(N) 

<  o  < 

N-l 


N-l 


(3.64b) 


and 


”,  <  '?> 


,  (N) 

<  \  '  '  <  \J 

N-2  N-l 


(3.64c) 


These  three  interlacing  conditions  are  related  to  each  other.  In  fact, 
(3.64c)  implies  (3.64a)  and  (3.64b).  One  can  see  this  by  writing  the 
quadratic  identity  as  follows: 


I 

o 


(v  2)  0  (v  2)  =  Y  (v  2)  K  (v  2 
n  on  on  on 


(3.65) 


2 

Then  (3.64c),  together  with  (3.63),  implies  that  K0(vn  )  alternates 
in  sign,  i.e.  (3.64a)  holds.  In  addition,  (3.64a)  and  (3.61)  imply 
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(3.64b).  Thus,  we  only  need  to  focus  our  attention  on  (3.64c).  This 

interlacing  condition  cannot  be  deduced  from  (3.63).  We  present  in 

the  appendix  an  example  of  a  trio  of  spectra  satisfying  the  small 

interlacing,  which  give  rise  to  a  X-spectrum  which  violates  (3.64c). 

N-l  N-l  N-l 

Therefore,  not  all  interlacing  sequences  {w  »  (v  and  ^n'l 

are  bona  fide  spectra  of  an  ' 1-bean .  (This  is  to  be  contrasted  with 

the  situation  for  an  N-string) .  To  be  such,  these  sequences  must 
satisfy  the  additional  constraints 


(-1)"  I  (v  )  <  0  ,  n-l,  2,  ....  N-l 
o  n 


(3.66) 


or,  in  terns  of  the  data. 


N-l  2  N-l  2 

"T-'T"  u.  9  i-r  v 

(l-  -L-r  (i-  -Sj) 

i=l _  vi  isH  ui 


(-!)"  y  -iit - i - ifJL 

^  1  N-l  2  N-l 


- —  u  .  — — —  u  . 

j=l  (1-  -^)  (1- 

i-1  W1  i^j 


The  stripping  procedure.  Let  us  assume  that  we  have  determined  IQ(u)  ) 

and  K  (u>2),  and  that  X^  and  thus  found  satisfy  all  the  necessary 

o  n  n 

interlacing  conditions.  The  next  step  consists  in  writing  (3.52)  for 
i=l  as  follows: 


o  (w2)  e.U2) 

— — —  =  {,  _  -1 — _ 

K  (w2)  °  K  (u>2) 

o  o 


(3.67a) 
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(3.67b) 


Y.  (u2)  =  Y  Ou2)  +  l  9  (ui2)  +  l  9  (uj2)  , 

1  o  o  o  o  1 


I  (uj2)  I,(u)2) 

— -  =  f  —± - 

2,  1  „  ,  2,  ’ 


Yx(«  ) 


Y1(oi  ) 


D^U2)  ,  D.  (J) 

o  2,1 

5~  ~  +  - 7“  » 


Y1(«  ) 


Y1(u  ) 


^(oi2)  =  Kq(oj2)  +  f  1D1(uJ2)  -  m1a)2Io(u2) 


By  dividing  -©  (w*-)  by  KQ(w  ),  which  are  two  polynomials  of  the  same 

degree,  the  quotient  is  a  constant  4  and  remainder  a  polynomial  of 

2  2  2 
degree  N-2  in  w  ,  namely  0^(ui  ).  Knowing  4  and  ©^(w  ),  we  can  find 

2  2  2 
Y^(uj  ).  Dividing  next  IQ(w  )  by  -Y^(w  ),  once  again  two  polynomials 

2 

of  the  same  degree,  we  infer  the  quotient  f  and  the  remainder  I^(w  ). 

2  2 

Similarly,  by  dividing  (co  )  by  Y^(u  )  which  is  of  lower  degree,  the 
2  2 

quotient  is  m^io  and  the  remainder  is  D^(u>  ).  Finally,  we  can 
2 

evaluate  K^(co  )  and  start  the  cycle  over  again. 

We  postpone  a  discussion  of  this  formal  construction  procedure 
until  later  and  simply  remark  at  this  stage  that  by  setting  u>-0  in 
(3.67d),  (3.67e),  (3.67a)  and  (3.67b)  we  would  get: 


(3.67c) 


(3.67d) 


(3 . 67e) 
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D1(0)  =  1  , 


v°>  =  -Fr  +  fi  =  -pi  • 


0,(0)  =  F,-e  f(n;  =  y  (x.-i  )  f.  , 

1  loo  ~  1  O  1 


(3.68) 


where 


Y, (0)  =  -F  +  2t  F  -  l  2F(N)  , 

1  2  o  1  o  o 


l  (x.-l  f  .  » 

5:  i  o  i 


fCN)  =  I  f.  . 
0  1  1 


1^(0)  can  be  deduced  from  the  quadratic  identity: 


hm  •  l  h  l  (vV2  h  -  <f  ‘W  V: 


These  expressions  are  similar  to  those  for  Dq(0),  ....  Iq(0)  except  for 
the  fact  that  the  first  segment  as  well  as  the  first  clothespin  are 
missing.  Thus,  by  means  of  this  procedure  we  have  "stripped-off"  a 
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little  segment  of  the  beam.  This  is  reminiscent  of  other  inverse  problems 
such  as  in  seismic  prospecting  wherr>  this  stripping-off  is  usually 
carried  out  in  the  time  domain  (Berkhout  &  van  Wulfften  Palthe,  1979). 

The  Stieltjes  theorem.  The  procedure  we  have  outlined  is  formal  since 
it  does  not  necessarily  yield  values  of  and  f  which  are  positive 

and  hence  physically  meaningful.  Indeed,  numerical  experiments  carried 
out  with  spectra  satisfying  the  grand  interlacing  have  revealed  that 
conditions  (3.41)  and  (3.66)  are  not  sufficient  to  guarantee  the 
positivity  of  the  physical  characteristics  of  the  N-beam.  In  contrast, 
the  interlacing  (3.2)  is  sufficient  to  guarantee  the  positivity  of  the 
density  of  the  corresponding  vibrating  string.  This  is  a  direct  result 
of  a  theorem  of  Stieltjes  which  we  can  state  thus: 

Part  1 .  If  0  <  6,  <  a  <  ...  <  ,  ,  then 

1  1  N— 1  N-l 


Part  2.  If  0  <  <  yl  <  ...  <  yn_2  <  ,  then 


I  I  (1-  ^“2  ) 

n=l  S 

n 


II  a-  ^ ) 

n=l  y 


I  I  (1-  ^2  ) 

,,  n=l  6  L 
2  n 

mu  +  - 


a-^  > 


n=l  Y 


(3.71) 


where 


m  >  0 


B .  <  6  ,  <  Y .  .  i=l,  2,  ...,  N-2  . 

i  i  i 


(3.72) 


We  can  apply  Part  1  of  Stieltjes  theorem  to  (3.67a)  to  show 


that  l  >0.  In  fact 
o 


0  F  (N)  v  2  ’ 

0  n=l  n 


(3.73) 


or,  equivalently 


(3.74) 
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Both  of  the  above  expressions  for  i  are  compatible  on  account  of  the 

2 

quadratic  identity.  Also,  if  we  denote  the  zeros  of  0^(w  )  by 
{v;2}?  2  ,  then  (3.70)  implies  that 


Vi  '  V  *  °i+l  *  i=1*  2.  ••••  N-2  • 


(3.75) 


Eliminating  l  from  (3.67b),  we  get 


2  2  2  2 

,  o  o'  (u )  0,  (W  ) 

Y.(u>  )  -  Y  (u  )  -  T—  +  ■■ ■  -r-  ■  , 

K  (u»  )  K  (u  ) 

o  o 


which,  because  of  the  quadratic  identity,  we  can  also  write  as 


?  I  (u>2)D  (to2)  Q. 2(w2) 

Y.  (oj  )  -  — - - 5  +  - r— 

K  (u»  )  K  (u)  ) 

o  o 


By  setting  cd=X^  in  the  above  equation,  we  deduce  that 


2  2 
2  Y<Xi  > 

Y  (X  )  -  -i - i—  ,  i-1,  2 . N-l 

1  1  K  (X.2) 

o  1 


Since  o^  and  X^  interlace,  K^(X^  )  changes  sign  as  i  goes  from  1  to  N-l, 

2 

Consequently  the  zeros  of  Y^(ui  )  and  X  interlace.  In  fact,  we 
can  show  that 


U1  A1  '  ’**  <  UN-2  <  xn-2 


(3.76) 
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This  is  the  condition  which  is  necessary  if  we  were  to  apply  Part  1 
of  Stieltjes  theorem  to  (3.67c).  As  a  consequence  it  follows  that 
f^  >  0  and  that 


X  .  < 

l 


X* 


i+1 


1=1,  2, 


N-3 


In  fact 


f(n)f  -f  2 

o  2  fl 


F  ’ 
2 


n 


n=l 


where 


F\ 


N 

=  l 
2 


(*i’V‘ 


A  more  useful  expression  for  f^  can  be  obtained  by  writing  (3.52e) 
thus : 


K1  " 


+  f,D  -  w  m,  I, 
1  o  11 


from  which  we  deduce  that 


f 


1 


N-l  2 

TT 

n»l  n 


(3.77) 


(3.78) 


(3.79) 
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Incidentally,  as  a  result  of  the  interlacing  between  w  and  o 

n  n 


(N) 


is  smaller  than  F 


(N) 


Hence  K^(0)  as  given  by  (3.68)  is  negative. 


(N) 

In  the  same  vein,  (3.73)  shows  that  IF  is  smaller  than  F,  and 

o  o  1 

hence  0^(0)  as  given  by  (3.68)  is  positive.  Finally,  in  this  same 


formula,  Y^(0)  can  be  seen  to  be  negative.  In  summary 


K:(0)  <  0  , 
01(O)  >  0  , 
Y1(0)  <  0  , 


i.e. ,  they  have  the  same  signs  as  Kq(0),  0q(O)  and  ^(0)  respectively. 
We  next  turn  to  the  equation  (3.67c)  for  and  proceed  to  eliminate 
f^  and  from  it.  In  other  words,  we  write  it  thus: 


K 


1 


Ii-1 

Ko  +  D1 


h  -D. 
o  1 


I 

o 


which  after  simplifications  brought  about  by  the  quadratic  identity 
for  i»l,  becomes 


-  D  I  +  Y  K  -  ©,  =  0  . 

o  o  1  o  1 


Setting  in  this  hybrid  quadratic  Identity,  we  see  that 


2  9lV> 

w2)  -  -  7777 

o  i 
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Making  use  of  (3.76)  we  deduce  that 


i  <  U  i  < 


“i+l 


1=1,  2,  N-2  . 


We  are  now  In  a  position  to  apply  Part  2  of  Stieltjes  theorem  to 
(3.67d).  As  a  result  m  >  0  and 


co .  <  to*  .  <  u\ 
ill 


i-1,  2, 


N-2  . 


In  fact 


m 


1 


N-2 


n=l 

N-l 


n=  1 


For  the  record,  we  also  write  an  alternative  formula  for  m  ,  nan 


(N) 


N-l 


ml  =  „ (N) 


p '  'r  _p 

o  '2  ‘  1 


I  ° 


(N) 1 


n=l 


We  can  derive  two  more  interlacings  by  setting  <d  equal  to  ar 


in  the  quadratic  identity.  These  are 


and 


<  X.’  <  o’ 

x  i+l 


1=1,  2,  ....  N-2 


(3.87) 


In  summary,  starting  with  the  polynomials  Io(u“),  ...,  DQ(u)  ) 
whose  zeros  satisty  the  grand  interlacing  (3.40)  we  can  deduce 

(i)  jIq ,  f  and  which  are  positive, 

2  2 

(ii)  new  polynomials  I^(ui~),  ...,  D^(w“)  of  lower  degree;  these 
polynomials  have  the  same  structure  as  the  original  ones  at  u=0.  How¬ 
ever,  their  zeros  satisfy  only  four  interlacing  conditions.  These 
conditions  are 


<  X, 


;  *  <-  .. » 

N-3  ^N-2 


<  X, 


Sl-3  <  °N-2 


(3.88) 


“N-2  < 


N-2 


<  a. 


^-2 


N-2 


Thus,  we  do  not  close  the  circle:  other  conditions  must  be  placed 
on  the  data  to  insure  that  the  eigenvalues  of  the  "stripped"  beam 
satisfy  the  grand  interlacing. 
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....  <y  V? 


!  in  i  am . 


Another  method  for  reconstructing  the  N-beam.  We  return  to  the 


continuous  beam  and  to  equations  (3.27’).  After  stripping-off  the 

portion  (o,x)  from  the  original  beam,  its  nth  eigenfrequency  in  the 

clamped/clamped  configuration  will  be  X  (x)  such  that 

n 


I(x,An"(x))  =  0  . 


(3.89) 


Varying  x  a  little,  we  see  that 


2  ?  dX 

I’(x,A  "(x))  +  i(x,A  (x))  ~n~ 
n  n  dx 


(3.90) 


2 

where  a  dot  represents  a  differentiation  with  respect  to  ui  .  But 
according  to  (3. 27 ’a) 


l(x. 


\>> '  - 


r2(L) 


X' 


(t-x)2 

r 


dt 


X' 


k^n 


(1- 


X  2(x) 

n _ 

Xk2(x) 


(3.91) 


and  I’  can  be  obtained  from  the  auxiliary  equations  (3.25). 
Consequently 
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dA 

n 

dx 


V(x) 

r  (>:) 


L 

dt 
J  r 


J  r 


X 

rL  2 

L  2 

(t'X)  dt-( 

'f  —  dt) 

r  \ 

k  J  r  / 

x  x 


X 


I  ! 

k^n 


1- 


u  (x) 

\2<*> 


v2<»> 

un2(x) 


1) 


Note  that  each  term  in  the  right  hand  side  is  positive.  Hence  X  (x) 

n 

an  increasing  function  of  x.  The  reader  may  recall  that  we  gave  a 
different  proof  of  this  result  earlier. 

In  the  same  manner,  we  can  show  that 


du 


.  2u  (x) 
dx  n 


dv 


t-x 


dt 


1- 


U  2(x) 
n 


vk2(x) 


dt  k#n  L~ 


U  2(x) 
n 


uk2(x) 


vn2(x) 


-J —  *  V  '■(x) 
dx  n 


dt 

r 


«k2(x) 


un2(x) 


v  2(x) 
n 


-  1 


k#n  1- 


v  2  (x) 

(  - - 1)  , 


°tT<x) 


Vk2(x) 
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n  v  2 ,  v  (t-x) 

— —  51  o(x)w  (x)  •  -  dt 

dx  n  r 


2  , 
w  (x) 

-p— i--V" 

»ic  <»> 

,  ^n2(x) 

k*n  "  u)k2(x) 


CO  2(X) 

(1-  -  ) 

un2(x) 


(3.96) 


Our  previous  discretization  for  the  N-beam  has  taught  us  (i)  that 

and  cannot  be  retrieved  and  (ii)  that  there  are  N-2  eigenvalues  for 
N 

the  clamped/clamped  configuration  as  opposed  to  N-l  for  the  others. 
Keeping  in  mind  these  lessons,  let  us  discretize  the  above  formulas 
thus: 
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CHICASO  UNIV  IL  DEPT  OF  GEOPHYSICAL  SCIENCES  F/6  20/11 

INVERSE  PROBLEM  FOR  THE  VIBRATING  BEAM  IN  THE  FREE/CLAMPED  CONF— ETCUJ) 
1979  V  BARCILON  N0001R-76-C-003R 

NL 


UNCLASSIFIED 


i 

\ 


If 


OJ 


,  2 
n 


oj  +  m,iu  F„ 
n  1  n  2 


N-l 


k=l 


N-l 


k^n 


(1-  ~^2  > 


(3.96’) 


The  reconstruction  of  the  N-beam  proceeds  as  follows: 

2 

(i)  by  means  of  the  interpolation  formula  (3.59)  for  Xq (uj  )  we 

infer  F  and  (X  ^  ; 
o  n  1  ’ 

2 

(ii)  from  the  quadratic  identity  we  get  Kq(u  )  and  hence 

,  ,N-1 

{an}l  ; 

(iii)  we  compute  f^  and  by  means  of  (3.74),  (3.79)  and 
(3.85); 

(iv)  we  compute  X*  ,  p’  ,  v’  ,  o’  and  w’  from  (3 .92’ )- (3 .96’ ) . 

n  n  n  n  n 

This  information  enables  to  continue  the  stripping  process. 

The  above  procedure  has  many  advantages  over  the  previous  one. 

It  is  attractive  from  the  numerical  point  of  view  since  long  divisions 

as  well  as  locations  of  zeros  are  avoided.  It  is  pleasing  from  the 

esthetical  point  of  view  since  it  couples  the  given  eigenvalues  to 

the  unknown  structure  without  the  intermediary  of  the  eigenvalue 

problem.  Finally,  it  might  prove  helpful  in  the  search  for  the 

missing  conditions  satisfied  by  bona  fide  data.  Indeed,  (3-92)— (3.96) 

indicate  that  the  relative  distribution  of  eigenvalues  determines  the 

slope  of  the  curves  w»X  (x) ,  ....  w-w  (x)  at  x=*0.  What  should  the  values 

n  n 

of  these  slopes  be  in  order  for  the  curves  not  to  intersect  each  other? 
The  answer  to  this  question  would  provide  the  missing  conditions  we  are 
looking  for. 
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'A 


i 


i 


The  limit  N-*».  We  would  like  to  examine  whether  or  not  the  sequences 
{p(N)(x)}  and  (r^(x)}  converge,  where 


(N)  ,  .  Nr1  (N)  x,  (N)  v 

p  (x)  =  l  m.  5 (x-x  ) 

i=l  1 


and 


(N)  ,  .  “  i 

r  (x)  i*l 


l  f.<S>  «,,«). 


We  shall  assume  that  the  beam  has  a  finite  length,  i.e. 


Un  (N) 
L  =  N-k°  xN 


or,  equivalently. 


where 


Un  N"1 

L  =  .Un  l 


(N) 


N-h»  ^  i 
i=*0 


*  (N)  -  £ 
o  Fx 


N-l 


2  ’ 


n*l  n 


(N) 


N-2 


V 

V 


n“l 


,•  2 
n 

.  2  ' 


etc. 


(3.97) 


(3.97’) 
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If  the  data  are  such  that  (3.97)  holds,  then  we  can  define 


F(x) 


x 


and 


M(x)  = 


(L-t)  p  (t)dt  . 


Clearly,  the  N-th  approximation  of  these  functions  are 


:(N) 


(x)  =  I  f£ 


(N) 


x^>x 


and 


.< (N)  ,  r  ,  (N)  (N) . 2  (N) 

M  (x)  *  i  -xi  )  mi 

X  .  >x 
1— 


From  their  definitions,  F^(x)  and  M^(x)  are  non- increasing  functions 
of  x.  In  fact 


F<N\x)  <  F  (N> 
—  o 


(N) 

We  can  express  Fq  in  terms  of  the  data  by  writing  (3.59)  as  follows 


,  (N)  _1_ 

o  "  F„ 


I  (0) 

o 


(3.98) 


(3.99) 


(3.98’) 


(3.99’) 


(3.100) 
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or  better  still 


N-l 

n 

2 

u 

1-  -^l 

2 

2  "> 

F/  F.z  N-l 

1  +  1  s 

k=l 

V1 

F,  F,  J-  N-l 

2  2  -1  pj 

2 

1-^ 
x  2 

N-l  | 

TT 

2 

l-!n_ 

2 

k=l 

-k 

k^n 

Therefore,  if 


(N) 

Then  the  limit  of  Fq  as  N  tends  to  infinity  exists,  i.e. 


F  -  lim  F  (N) 
o  N-*»  o 

(N) 

and  this  limit  provides  an  upper  bound  for  F  (x) ,  i.e. 

(N) 

F  (*)  1  F0  for  0  £  x  1  L  and  for  every  N. 


(3.101) 


(3.102) 


(3.103) 


(3.104) 


(N) 

Similarly,  the  functions  {M  (x) }  are  also  bounded  above.  We  can  see 

this  as  follows.  Let  us  substitute  the  series  representations  for 
2  2 

Y(x,u>  )  and  D(x,u>  ),  namely 


3.45 


L 


Y(x,u)2) 


-r2(L)  [  J  ilf 


dt  +  £  (-l)n  y  (x)o 


(x,u2)  =  r2(L)  [  1  +  l  (-1)"  dn(x)o>2n  ] 
*■  n=l 


in  the  auxiliary  equation  (3.25f).  Then,  the  linear  terms  in  m  yield 


L  L 

(t-t')2 

d1(x)  -  p(t')df  dt 


Confronting  this  result  with  the  product  representation  of  D(x,u>  ) 
given  in  (3.27’e),  we  conclude  that 


co 

y  -1_ 


o(t')dt' 


f  (t-t'): 

j  r  (t) 


or,  after  interchanging  the  order  of  integration 


(3.105) 


c 


00 


L 

( 


J  r (t) 

0  0 


2 

(t-x)  p(x)dx. 


(3.106) 


— being  integrable  for  t  in  (0,L)  we  conclude  that 

t 

2 

(t-x)  p(x)dx  <  °°  . 

0 

As  a  result  A/(0)  exists  and 

(x)  £  iV(0)  for  0  x  £  L  and  for  every  N.  (3.107) 

(N) 

Therefore,  the  non- increasing  sequences  of  function  [F  (x) }  and 
(N) 

{ M  (x)}  are  bounded  above.  Hence,  according  to  Helly's  selection  theorem 

theoren  (see  e.g.  Natanson  1955,  p.  220),  we  can  extract  from  these 
sequences  two  subsequences  which  converge  for  all  x  in  (0,L)  to  two 
non- increasing  functions.  In  view  of  the  uniqueness  theorem,  these 
limit  functions  must  be  the  functions  F(x)  and  M(x)  associated  with 
the  non-pathological  beam.  Finally,  p(x)  and  r(x)  are  obtained  after 
differentiation. 


3.47 


CONCLUDING  REMARKS 


We  saw  that  the  information  contained  in  the  impulse  response  is 
equivalent  to  the  three  symphathetic  spectra  {wn},  {v^}  and  { wn } *  and  the 
two  moments  of  the  flaccidity  and  ?£.  We  also  saw  that  the  information 
contained  in  the  impulse  response  does  not,  strictly  speaking,  determine  a 

single  beam.  Rather,  it  determines  a  class  of  beams:  all  the  beams  in  this 

class  have  the  same  oscillating  portion,  but  they  differ  over  that  rear  portion 
of  their  lengths  which  is  stationary.  This  stationary  portion  is  made  up  of 

a  weightless  rod  of  infinite  rigidity.  I  was  surprised  by  the  fact  that  the 

massive  wall  in  which  the  beams  are  embedded  is  not  part  of  the  ambiguity. 

The  uniqueness  proof  has  the  drawback  of  requiring  that  p(x)  and  r(x) 
be  differentiable  functions.  This  is  doubly  regrettable:  firstly,  because 
the  degree  of  smoothness  of  p(x)  and  r(x)  is  not  easy  to  infer  from  the  data 
and  secondly,  because  discontinuities  are  likely  to  occur  in  the  more  difficult 
geophysical  problem  for  which  the  beam  is  to  serve  as  a  guide.  Ideally,  these 
smoothness  requirements  ought  to  be  relaxed.  To  that  effect,  it  would  be 
desirable  to  derive  the  uniqueness  result  from  the  construction,  i.e.  to  show 
that  the  approximations  form  a  Cauchy  sequence. 

Our  discussion  of  the  inverse  problem  was  carried  out  entirely  in  the 
frequency  domain.  Yet,  the  stripping  procedure,  which  unraveled  progressively 
the  structure  of  the  beam,  is  a  reminder  of  the  presence  of  a  wave  propagating 
along  the  beam.  In  spite  of  this  time  domain  intrusion,  I  believe  that  the 
frequency  approach  is  best  suited  for  this  inverse  problem  at  hand.  Indeed, 
the  continued  fractions  reliance  on  the  w-dependence  shows  very  clearly  the 
advantage  of  working  in  the  frequency  domain.  In  addition,  the  theories  of 
entire  functions  and  of  oscillating  kernels  provide  some  very  powerful  tools 
for  the  investigation  of  the  problem  in  the  frequency  domain. 


4.1 


The  actual  construction  of  the  solution  required  a  discretization  of 
the  original  problem.  This  discretization  is  an  essential  step  which  cannot 
be  by-passed.  Indeed,  the  global  condition  (3.102)  is  a  direct  result  of  this 
apparent  detour.  The  same  remark  holds  for  the  vibrating  string:  the  global 
condition  (3.3)  is  a  consequence  of  a  discretization  followed  by  the  limiting 
process  N-^. 

The  existence  of  the  solution  to  the  inverse  problem  is  tied  to  very 

stringent  interlacing  conditions  on  {at  },  {v  }  and  {p  }.  These  conditions 

n  n  n 

are  such  as  to  guarantee  that  the  spectra  {w  (x)},  ....  {X  (x) }  for  the 

n  n 

truncated  beam  spanning  the  interval  (x,L) ,  satisfy  the  grand  interlacing. 
Rather  than  having  to  check  at  each  step  of  the  construction  whether  the 
grand  interlacing  is  satisfied,  it  would  be  desirable  to  have  explicit  ways 
of  testing,  from  the  outset,  whether  the  data  are  bona  fide. 

In  the  investigation  of  the  limiting  process  N-»®,  the  length  L  of  the 
beam  was  assumed  to  be  finite.  This  convenient  assumption  is  not  essential. 
Very  minor  modifications  are  needed  to  handle  the  infinite  case  since  the 
Helly  selection  theorem  can. still  be  used. 

Finally,  our  analysis  dealt  exclusively  with  the  beam  in  the  free/ 
clamped  configuration.  For  other  vibrating  configurations,  the  details  of 
the  results  would  be  modified:  for  instance,  the  impulse  responses  would  be 
equivalent  to  other  trio  of  sympathetic  spectra  and  the  gross  constants  would 
not  necessarily  be  related  to  moments  of  the  flaccidity.  However,  the  approach 
to  the  questions  of  uniqueness,  existence  and  construction,  which  was  used  in 
the  present  paper,  ought  to  be  applicable  for  general  boundary  conditions. 
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FIGURE  CAPTIONS 


Figure  1: 


Figure  2: 


Figure  3: 


Figure  4: 


Sketch  of  the  Green's  function  N(x,t)  for  the  non  self-adjoint/ 
clamped  configuration. 

The  curvature  of  N(x, t)  changes  at  x  =  xQ(t). 

2  ? 

Sketches  of  the  auxiliary  function  I(x,a)  )  for  a  value  of  oj” 

2  2 
lying  between  A^  and  A 2  • 

The  number  of  local  extrema  depends  on  the  relation  between  the 
A-  and  u-  spectra. 

2  2  2 
The  to  -x  plane.  The  solid  curves  A^  (x)  and  A?  (x)  represent  the 

first  and  second  eigenvalues  of  the  truncated  beam  spanning  the 

interval  (x,L)  in  the  clamped/c lamped  configuration.  The  dotted 

curves  represent  the  corresponding  eigenvalues  for  the  supported/ 

clamped  configuration.  By  considering  I(x,A2  )  we  can  show  that 

this  arrangement  of  eigenvalues  is  not  permissible. 

The  discrete  N-beara.  Its  structure  consists  of  N  clothespins 
located  at  x^,  X£,  ...»  x^.  These  clothespins  have  masses  {m^} 
and  flaccidity  {f_^}  and  are  connected  to  each  other  by  weightless, 
infinitely  rigid  rods  of  length  { P. ^ } . 
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appendices 


Alternative  discretization.  We  still  insist  that  the  mass  of  the 
beam  be  concentrated  in  a  few  points,  namely  that 


N 

ci  (x)  *  l  m  5(x-x  ) . 

i=l 


However,  we  now  allow  the  connecting  segments  to  have  a  finite 
rigidity: 


1 

r(x) 


for  x.  ,  <  x  <  x.  . 
l-l  —  —  1 


Substituting  these  expressions  for  o(x)  and  r(x)  in  the  auxiliary 
equations  (3.25  ),  we  deduce  that  over  the  interval  (x^  x^) 

D(x ,u>2)  =  Di_1("2)  , 

2  °t-l(w2)  2 

K(x,uj  )  =>  -  •  (x-x.  )  +  K.  (w  )  , 

i-1  1  1_J- 

7  D1-1(w2)  2  7 

Q(x,wZ)  *  -  •  (x-xi_1>‘i  +  )  •  (x-x^) 

+  * 


L 


A.  1 


These  are  the  discrete  auxiliary  equations.  It  should  be  noted  that 
they  satisfy  the  following  identity: 


-XiDi  +  YiKi  - 


-ri-lDi-l  +  Wi-l 


This  identity  leads  to  the  quadratic  identity  when  coupled  with  the 
boundary  conditions. 


Example  of  a  violation  of  the  grand  interlacing.  Let  us  consider  a 
4-beam  whose  impulse  response  satisfies  the  small  interlacing.  More 
specifically,  let 


u^2  *  1.00,  Vj2  =  2.00,  u22  =  3.00  , 

u>22  =*  3.01,  v2“  =  4.00,  u 2*"  =  3.00  , 


6.00, 


7.00, 


=  8.00 


and 


Fx  »  1.00 


F2  =  1.33  . 


2 

Following  the  method  outlined  in  §3,  we  can  construct  Iq(u)  )  which  is 

2  2 

a  polynomial  of  degree  2  in  u>  .  The  zeros  of  IQ(w  )  turn  out  to  be 


4.94 


Note  chat 


.  2  2 
A1  >  v2 


,  2  .  2 
X2  >  v3 


which  violates  the  grand  interlacing.  Were  we  to  pursue  the  construction 
of  this  4-beam,  we  would  find  that  some  of  the  i's,  m's  and  f’s  are 


negative. 


